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These introductory notes are about functional renormalization group equations and some of their 
applications. It is emphasised that the applicability of this method extends well beyond critical 
systems, it actually provides us a general purpose algorithm to solve strongly coupled quantum 
field theories. The renormalization group equation of F. Wegner and A. Houghton is shown to 
resum the loop-expansion. Another version, due to J. Polchinski, is obtained by the method of 
^^f^ ' collective coordinates and can be used for the resummation of the perturbation series. The genuinely 

' non-perturbative evolution equation is obtained by a manner reminiscent of the Schwinger-Dyson 

, equations. Two variants of this scheme are presented where the scale which determines the order 

^SJ ' of the successive elimination of the modes is extracted from external and internal spaces. The 

—^ \ renormalization of composite operators is discussed briefly as an alternative way to arrive at the 

^ ■ renormalization group equation. The scaling laws and fixed points are considered from local and 

h— ^ ' global points of view. Instability induced renormalization and new scaling laws are shown to occur 

, in the symmetry broken phase of the scalar theory. The flattening of the effective potential of a 

I ■ compact variable is demonstrated in case of the sine-Gordon model. Finally, a manifestly gauge 

invariant evolution equation is given for QED. 
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I. INTRODUCTION 



The origin of renormalization goes back to hydrodynamics, the definition of the mass and other dynamical charac- 
teristics of bodies immersed into fluids. The more systematic elaboration of this concept is based on the (semi) group 
property of changing the observational scale in a theory. The renormalization group (RG) method has already been 
used with a number of diflferent goals. Some of the more important directions are (i) to remove U.V. divergences [1], 
(ii) to describe the scale dependence of physical parameters and to classify the parameters of a theory around a critical 
point according to their impact on the dynamics [2] , (iii) to express the highly singular product of local field variables 
[3] , (iv) to resum the perturbation expansion in Quantum Field Theory [4] and in the case of differential equations [5] 
and finally (v) to solve strongly coupled theories. Furthermore, the RG method offers a systematic way of organizing 
our understanding of a complicated dynamics by successively eliminating degrees of freedom. The distinctive feature 
of the method is that it retains the influence of the eliminated subsystem on the rest. 

The main subject of these lectures is application (v) in the framework of the functional formalism, the search of a 
general purpose algorithm to handle non-perturbative and not necessarily critical theories. 

The strategy of the RG scheme for the resummation of the perturbation expansion is to introduce and to evolve 
effective vertices instead of dealing with higher order corrections. The result is an iterative procedure where the 
contributions to the effective vertices computed at a given step in the leading order of the perturbation expansion 
and inserted in a graph at the next step reproduce the higher order corrections. In the traditional implementation of 
the RG procedure one follows the evolution of few coupling constants only. Our concern will be an improvement on 
this method in order to follow the; evolution of a large number of coupling constants. This is realized by means of the 
functional formalism where the evolution is constructed for the generator function for all coupling constants. 

Before starting let us review the conventional multiplicative RG schemes as used in Quantum Field Theory. This 
scheme is based on the relation 



G^') (pi , • • • , Pn\ 9R, m) = Z-^/^G^^^ (p,, . . . A) + O ( ^ 




(1) 



for the Green functions, where A and /x denote the U.V. cut-off and the substraction (observational) scales, respectively. 
Z depends on either gji or gs and the ratio A//i (the mass is treated as one of the parameters g). The ignored terms 
stand for non-universal, cut-off dependent interactions. 

Bare RG equation: The renormalized theory is independent of the choice of the cut-off. 
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where the renormalized coupling constants and the substraction scale are kept fixed in computing the derivatives. 
The equation involving the renormalized quantities is not too useful. It is the multiplicative renormalization scheme 
which leads to the second equation, expressing the possibility of compensating the change of the substraction scale by 
the modification of the coupling constants gsiA) arising from the solution of the differential equation by the method 
of characteristics. Notice that the compensation is possible only if the list of the coupling constants {93} includes all 
relevant coupling constants of the appropriate scaling regime, around the U.V. fixed point. 

Renormalized RG equation: The bare theory is independent of the choice of the substraction scale. 



• • ,Pn;5i?,M) = Z "^'^i9R,IJ')GR\pi,--- ,Pn;9R,fJ') 



0, 



(3) 



where the bare coupling constants and the cut-off are kept fixed. 

Callan-Symanzik equation: The change of the mass in the propagator is governed by the expression 
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(4) 



which can be used to find out the dependence on the renormalized mass when the cut-off and all bare parameters 
except the bare mass are kept fixed. The resulting evolution equation [6] is similar to the two previous RG equations 
except that the derivative is with respect to the bare mass instead of the cut-off or the substraction scale and the 
right hand side is non-vanishing. This latter feature indicates that contrary to the first two schemes the evolution 
9R{m?) in the Callan-Symanzik equation is not a renormalized trajectory, it connects theories with different mass. 

The serious limitation of these equations is that they are asymptotic, i.e., are applicable in the regime iJ.^,p^ « A^ 
only. In fact, the omission of the non- universal terms in the multiplicative renormalization scheme (1) requires that 
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wc stay away from the cut-off. In models witli IR instability (spontaneous symmetry breaking, dynamical mass 
generation,...) another limitation arises, Tri^y^ << iJ.^,p^, to ensure that we stay in the U.V. scaling regime. This is 
because the IR scaling regime may have relevant operators which are irrelevant on the U.V. side [7]. In more realistic 
models with several scaling regimes which have different relevant operator sets we need non-asymptotic methods 
which can interpolate between the different scaling regimes. This is achieved by the functional extensions of the RG 
equation based on the infinitesimal blocking steps. 



II. FUNCTIONAL RG EQUATIONS 

There are different avenues to arrive at a functional RG equation. The simplest is to follow Wilson-Kadanoff 
blocking in continuous space-time. It leads to a functional differential equation describing the cut-off dc;pcndcncc of 
the bare action. Better approximation schemes can be worked out when the modes are suppressed smoothly as the 
cut-off changes. The parameters of the bare action which are followed by these RG equations are related to observables 
qualitatively only. It is more advantageous to transform the RG equation for the effective action whose parameters 
have direct relation with observables. 

For the sake of simplicity we consider an Euclidean theory in d dimensions for a scalar field (t)x governed by the 
action S's [(/)]. An 0(d) invariant U.V. cut-off A is introduced in the momentum space by requiring (/>p = for \p\ > A 
to render the generator functional 

ei^y^ = J D[<f>]e-i^''Mr'l> (5) 

finite. We shall use the following notation conventions, stated below: In order to render the functional manipulations 
well defined we always assume discrete spectrum, i.e., the presence of U.V. and IR regulators, say a lattice spacing 
a and system volume V — a'^N'^. The space-time integrals are J d'^x = n.'^J^x = Jx> Ix f^S^ ~ -f ' -9' / d'^'p/ {'^''^Y = 
y-^ = /p> fp = Ix e"*^^/a:, and fx = e'P'^ fp. The Dirac-deltas 5x,y = a^'^d^y and (5p,g = Vd^^ are expressed in 
terms of the Kronecker-deltas 5^,, and 5^„. 



A. Resumming the loop expansion 

We start with the simplest form of infinitesimal blocking-step RG equations which is the functional extension of 
the bare RG scheme mentioned above [8]. 



1. Blocking in continuous space-time 

We shall denote the moving U.V. cut-off by k. Its lowering k k — Ak leads to the blocking transformation of 
the action which preserves the generator functional (5). Due to the presence of the source this blocking introduces an 
explicit source dependence in the action, 

J DkW-^^'^^'t'Mr'^ = j Dk-i^kW'^^"-^"^''''^^^^^"^ , (6) 

where Df.[4>\ stands for the integration measure over the functional space Tk consisting of functions whose Fourier 
transform is non- vanishing for \p\ < k. In order to avoid this complication one usually assumes j G J^k-Ak- In this case 
the blocking becomes a mapping Sk[(f)] — > Sk-Ak[(l>] and it is enough impose the invariance of the partition function 

e-*Sfc-AfeM = /"£)[0]e-*S'=['?^+<^l, (7) 



where (/) e J^k-Ak and (j) & ^k\^k-Ak- 

The evaluation of the path integral by means of the loop-expansion gives immediately the functional RG equation, 

Sk-Ak[^] =Sk[^ + U + ^tr In '^'^^[^ + ^-^] + o in') , (8) 
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where 4>ci denotes the saddle point. The trace is over the functional space 4> S J^k\-^k-Ak and O represents the 
higher loop contributions. We write this equation as 



Sk[^] - Sk-Ak[^] = Sk[ct>] -Sk[^ + ici] - Jtr In ^ ^"J^t ^''^ + O {fi') . (9) 



tree level 



loop contributions 



When can one safely assume that the saddle point is trivial, 4>ci = 0? Let us suppose that (j)x is weakly and slowly 
varying, ie = ^ + r]x where r]x ^ and the characteristic momentum of r] is small with respect to the cut-off k. 
Then the fluctuating component rj appears as an external source breaking translation invariance. The saddle point, 
4> G ^k\^k-Ak, being inhomogeneous exists as 77 ^ only if the external space-time symmetry is broken dynamically. 

But a Afc-dependent small saddle point, ^ci = 0(AA;"), n > may occur without breaking external symmetries. 
In this case we expand in ^ci and find the Wegner-Houghton (WH) equation [8] 

^.[<^]-5.-A.M=2^-^^J (10) 
_^trln4^ + 0(?i')+0(Afc2) 



where all functional derivatives are taken at = 0. The saddle point can be omitted in the argument of the logarithmic 
function because the trace brings a factor of Afc as we shall see below. The discussion of the tree-level renormalization 
when (pel is non- vanishing and finite as Afc ^ is deferred to section IIIC. 

One can understand the loop contribution better by splitting the action into the sum of the quadratic part and the 
rest, S = ■ Gq ^ ■ (/) + Si and expanding in St, 

S(I>S<I) ^1 " V S(j}S(l) J ^ ' 

We recovered the sum of one- loop graphs. The loop corrections close on all possible pair of ^ legs with the propagator 
Go- The tree- level piece describes the feedback of the change of the cut-off on the dynamics of the classical background 
field 0. This is classical contribution because the cut-off controls the configuration space, the number of degrees of 
freedom. Some of the graphs contributing to the right hand side are shown in Fig. 1. 



(a) 



(b) 

FIG. 1. Graphs contributing to the blocking, (a): tree-level, (b): one-loop corrections. The dashed line stands for a particle 

of momentum |p| = fc and the solid lines represent <p. 

The point of central importance is that AA; serves here as a new small parameter to suppress the higher-loop 
contributions. As an example consider the simplest case, = <i>, (/)c; = with the ansatz 



SkW\ = f 

J X 

for the action with Z{(j)) = 1, 



\Zk{cl>x){d^cl>xf + Uk{(l>x) 



(12) 



Uk-Ak{^) = Uk{^) + \ j \n\p^ + U'^m+0{{hAkf). (13) 

^ Jk-Ak<\p\<k 
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Each loop-integral is over the shell k — Ak < \p\ < k in momentum space. So long the propagator is non-singular in 
the integration domain the n-loop integrals will be proportional to the n-th power of the integration volume, giving 
a dimensionless small suppression parameter « (AA://c)",the proportion of the modes eliminated with those left. The 
question of the singularity will be considered in section IIIB2. The higher loop contributions are suppressed in the 
infinitesimal blocking step limit and the one-loop evolution equation turns out to be an exact functional equation. 
The limit Ak ^ is safe for the loop contributions but more care is needed when the saddle point is non-trivial 
because the tree-level contribution might be singular in this limit. 

Notice that the convergence of the loop expansion was assumed in the argument above. Another word of caution 
is in order about Eq. (9). Not only the trace but the logarithm function itself is considered within the subspace 
4> € ^k\^k-Ak, a rather formal observation which becomes essential later. 

It is sometimes easy to arrive at exact, but unmanageable equations. The problem is rather to find an approximation, 
an acceptable compromise between precision and simplicity. By assuming the absence of non-local interactions and 
the homogeneity of the vacuum the blocked action is usually expanded in the gradient of the field, leading to the 
ansatz (12) for the action where contributions O (9^) are neglected. This step is the key to the success of the RG 
method. The point is that the blocked action can be identified by evaluating the integral (7) for different background 
configurations. The determination of the blocked action is an over-constrained problem because there are 'more' 
background field configurations than coupling constants in a given ansatz. The tacit assumption, that the effective 
action is well defined, i.e., independent of the set of configurations used to read off its coupling constants, gives the 
real power of the RG method: it enables us to make predictions. 



2. Local potential approximation 

In order to relate the WH equation to conventional perturbation expansion we shall use the ansatz (12) with Z = \, 
in the local potential approximation. In order to determine the only non-trivial piece of this action the potential 
Uk{(f>) it is sufficient to consider a homogeneous field 4i{x) = $ in (9). The saddle point is vanishing as mentioned 
above, (pel = 0, so long as the external space-time symmetries are unbroken. The blocking equation simplified in this 
manner yields Eq. (13), the projection of the WH equation into the restricted space of actions, 

^''(^^ = -^Mk'+ui'\n (14) 

where the dot stands for kdk and flj, = 27r''/^/r(d/2) denotes the solid angle. 



FIG. 2. The first four one-loop graphs contributing to the WH equation in the local potential approximation when the 
potential U{(t>) is truncated to the terms cjy^ and 0*. The dashed line corresponds a particle of momentum \p\ = k. 

The expansion of the logarithmic function leads to the series 

where — U^^\0), the sum of the one- loop graphs whose external legs have vanishing momentum and the internal 
lines carry momentum \p\ = k, cf Fig. 2, the leading order contributions to the renormalization of the coupling 
constants g„ in the loop expansion. 

The partial resummation of the perturbation series performed by the WH equation can easily be seen by comparing 
the solution of Eq. (14) with those obtained in the independent mode approximation, where the fc-dependence is 
ignored in the right hand side of the equation, 
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Ukm = Uxm-^ [ ln[p2 + [/f ($)]. (16) 

^ Jfe<|p|<A 

One recovers the one-loop effective potential at the IR fixed point, K;ff($) = J7fe=o(^)- This is not by accident and 
remains valid for the complete solution, as well, because e"^^""^*' is the distribution of the homogeneous mode and 
the blocking (7) at fc = retains the homogeneous mode only, S'fe=o[*^'] = VUk=o{^)- Wc shall see in section IIIC 
that this simple argument is valid in the absence of large amplitude inhomogeneous instabilities only. 
Let us consider as an example the differential equation 

ifc = /(ajfc), xx=xb, (17) 

where f{x) is weakly varying analytic function given in its expanded form around the base point x = 0, 

(18) 

n 

The RG method allows us to perform the expansion at the current point, 

Xk-Ak =Xk- Akfixk) + O (Afc2) , (19) 

cf Eq. (8). The fc-dependence of the right hand side represents the accumulation of the information obtained during 
the integration. In a similar manner, the virtue of the functional RG scheme (9) is the expansion of the path integral 
around around the current action, instead of the Gaussian fixed point. 

A more detailed picture of the RG flow can be obtained by inspecting the beta functions. For this end we parametrize 
the potential as 

Uk{<i>) = J2'-^i^-^or, (20) 

n 

and write the beta functions corresponding to (po = {(j)) as 

/3„ = kdkgnik) = kdkd;Uk{^o) = d^kdkUki^o) = '^^^^ ln[e + djUkm (21) 



It is easy to see that /?„ is an n-th order polynomial 



(22) 



of the expression Gn = gn/{k^ + §2), eg 

P2 = G4 — G3, 

= - 3G3G4 + 2Gl 

Va = Gq- AG3G5 + UGjG^ - SGl - 6Gt, (23) 

= G7 - 5G3G6 + 20G^G5 - IOG4G5 - 6OGIG4 + 3OG3GI + 24Gl, 
Vg = Gs- 6G3G7 + SOG^Ge - SGrGr - 120G:^G5 + I2OG3G4G5 
-270G|G^ - lOG^ - IOG4G6 + 360G^G4 + 30Gl - 120G|, 

etc. The correspondence between these expressions and Feynman graphs contributing to the beta functions in the 
traditional renormalization scheme is indicated in Fig. 3. 





(a) 
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(b) 

2 

(c) 

FIG. 3. Graphs contributing to (3n, (a): n = 2, (b): n = 3 and (c): n = 4. 

We shall need later the beta functions corresponding to coupling constants made dimensionless by the running 
cut-off. By means of the parametrization </> = k^'^~'^^/^^, = 

bn]=rf + n(l-0, (24) 

U{(t>) = k'^U{4>) one finds 

/3„(0) = tnik) = k-^'-^0r,{ct>) - [9n]9n = "^^^"[Gs, ' " ' , G„+2] - [5n]5n, (25) 

where G„ = + 92)- 

An instructive way to read Eqs. (20) and (14) is that the potential J7(^) is the generator function for the coupling 
constants. The functional RG method is an economical book-keeping method for the computation of graphs and their 
symmetry factors by keeping track of their generator function(al). 

The higher loop contributions to the running coupling constants are generated by the integration of the differential 
equations kdkQn = The us consider the two-loop graph depicted in Fig. 4 and expand the subgraph in the square 
in the momenta of its external legs. The zeroth order contributions is the last graph of Fig. 3b. By going higher 
order in the gradient expansion one can, in principle, generate the full momentum dependence of the subgraph. In 
general, the n-th loop contributions appear in the integration after n step k ^ k + Afc and all loops are resummed in 
the limit AA; 0. 




FIG. 4. A two- loop contribution to §2- 




3. Gradient expansion 



The natural expansion for the long distance properties of systems with homogeneous ground state is the expansion 
in the inverse characteristic length of the fluctuations, in the gradient operator acting on the field variable. 
Functional derivatives of local functionals: We consider the following generalization of the O (9^) ansatz (12), 



S[ct>] 



(26) 



We shall use the notation = Z{(j)x), Ux = U{(j)x) for the coefficient functions and Z^"-* = d^Z{(j)x), J/i"^ = d^U{(j)x) 
for their derivatives. The first two functional derivatives of the action are 
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v— = ye'py— 



J x,y 



and 



V 



2 5^S 



fpuyq Jx,y,z 



Jpy+iqz ^ 



+ 5x,z {d^Sa:,yK~^dij,(px + d^(pa:K~^dij,Sx,y 
+6x,y6x,zUP^- 

We split the field into a homogeneous background and fluctuations, ?!>x = ^ + and find up to O (??^), 

+{Z + Z<-'^r],) {d^S,,yK-'d^r], + d^r,,K-'d^6,,^ 

J x^y I \ Jr ^ Jr,s / 

J r,s 

+ (^Z + ^ e'^-^rjr^ jp- se'^-^rts {Kj^ + R-') 
+1 [ 6p+r+s,oVrVs [t^^^^ - Z^^^ {r ■ sK-^ + p ■ rK-^ + p ■ sR-^) 

^ Jr.s '- 



and 



6^S 

(ppHq Jx,y,z 
1 



Jpy+iqz 



Sx,ySx,z (U^'^+U^%X + IU^'^VI) 



Z^^'^di_,r],K-'d^Vx 



+ (^Z + ^(1)77, + \z^^^if)j {d^5,,zR-'di,6,,y + d^6,,yR-^di,S,,,) 



(27) 



(28) 



(29) 



Jx,y,z I t/r,; 



+ t/(2) + 



'p+g,o 



C/(3)_V)(^p.g(;f-l+i^-l) 



Sp+q+r+sfiVrVs 



+q ■ s{K-^ + K-^) +p- r{K-^ + K;^) 



(30) 



where = = and i^^ = K{s^). 

WH equation in C (rf) : We set iiT = 1 in (26), use sharp cut-off, i.e. all momentum integral is for k — Ak < \p\ < k, 



S^S 

5<j)p5(j)q 



+ 9 / ^p+q+r+sfi'Hr'ns 
^ Jr,s 



[/(3) _2'(i)(p.g + g.r+p-r) 



(31) 



with 



— Sp+qfi 



U(^) + \z{p'+q^) 
= I Sp+<,+r,oVr C/'"' + Iz^'^ip'' + q^+r^ 

Bp,q = i ^ 5p+q+r+sfiVr1ls U^"'^ + \Z^'^\P' + 9^ + + s2) 



(32) 



As emphasised above, the logarithmic function on the right hand side of the WH equation is computed within the 
space of field configurations whose Fourier transform is non- vanishing for k — Ak < |p| < fc only. In order to keep this 
constrain in an explicit manner we introduce the projection operator V corresponding to this space and write the rest 
of the O (jy^) evolution equation as 



dkS[cf>] 



-— trlog[G-+A + B] 

^ ;trPlogG-i - ^tv[VG{A + B)] + J^t-,[VGAVGA]. 



2Ak 



4AA; 



(33) 



Since the propagator does not lead out from the function space, [P, G] = 0, it is sufficient to make the replacement 
G VG on the right hand side of the evolution equation. The term O {r}) drops owing to the vanishing of the 
homogeneous component of the fluctuation, r/p=o = and we flnd 



10 



2 / V-rVr 



2(27r) 

- / V-rVr 
J r 



2 del Flog[ZA;2 + C/(2)] 



r—ek.ek—r ~ 



[/(3) + izW(fc2 + (eA;-r)2+r2) 



+ 



j V-rVr 



2(Zfc2 + UC^)){Z{ek - r)2 + C/(2)) 

2(zfc2 + «7(2)) r 



(34) 



including an integration over the direction of k, the unit vector e. 
The identification of the terms O (r]^r°) results the equations 



U 



(35) 



c.f. Eq. (14). The terms O [rfr'^) give 



2(2- 

'T^r—ek.ek—r 



ZA;2 + C/(2) 



^^"^^ [zk-' + m^) (zfc2 + c/(2)) 




(36) 



Finally, in O (rj'^r^) we find 



(ZWfc2 + C/(3)) 



21 



ZA;2 + C/(2) 



(37) 



Let us make the replacement Vr-k.k-r ^ l + [7'r-fe,fe-r^l] in Eqs. (35)-(37) and call the contributions corresponding 
to 1 and Pk-r,k-r — 1 regular and irregular, respectively. The regular contribution can be obtained from the one-loop 
graphs. This has already been seen for Eq. (35) and the regular part of Eq. (37) is just the second derivative of Eq. 
(35). 

But there are problems with the irregular contributions which represent the cut-off in the multi-loop integrals. One 
obvious problem is the inconsistency between Eqs. (35) and (37). Another problem is that the irregular contributions 
include a truncated r-integration and as a result the left and right hand sides of the equations have different rj- 
dependence. This could be avoided by imposing r << Ak for the momentum of the fluctuations. But the price is 
unacceptable high since it would indicate that the radius of convergence for the gradient expansion is smaller then 
the infinitesimal AA;. 

Another problem the projection operator leads to is non-locality. It is known that the bare action is non-local at 
the scale of the U.V. cut-off, i.e. it contains higher order derivatives up to a finite order. In fact, the gradient 
represents couplings up to n lattice spacing distances when lattice regularization is used. But terms in the action 
which are non-polynomial in the gradient induce correlations at finite, i.e. cut-off independent distances and are not 
acceptable. Returning to the gradient expansion we note that the operator Vr-k,k-r restricts the integration domain 
for r in such a non-isotropical manner that not only but |r| appears, too. This latter corresponds to the non-local 
operator ^/^ in the gradient expansion. 

This problem, inherent in the sharp cut-off' procedure in the momcntiim space is rather general. The one-loop 
correction for the two-particle scattering amplitude in the (t>* model, represented by the fourth graph in Fig. 3c has 
a non-local contribution when the momenta in both internal lines are properly restricted by the U.V. cut-off. Such 
non-local effects, disregarded in the usual textbooks, might well be negligible in a renormalized theory because they 
are represented by non-renormalizable, irrelevant operators such as (p^^/Oip. Therefore one hopes that their effects 
become weak when the cut-off is removed. The more careful analysis is rather involved since the details of the cut-off 
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procedure within the momentum regime (1 — e)fc < p < (l + e)fc influences the dynamics on the distance scale h/ {2ek) 
which is a quantity of the form • oo. 

But it is not so easy to dispel the doubts. First, the irrelevance and unimportance of a coupling constant are 
quite different concepts. An irrelevant coupling constant approaches its I.R. fixed point value as we move in the 
I.R. direction. But its value at the fixed point may be strong, indicating that the coupUng constant in question is 
important. It is the scale dependence only what becomes unimportant for an irrelevant coupling constant and not its 
presence. The higher order vertices induced at low energies in a strongly coupled self interacting scalar field theory 
in < 4 may serve as an example of this difference. Second, the global view of the renormalization group, outlined in 
Section III B suggest that a coupling constant which is irrelevant at the U. V. fixed point docs not necessarily remain 
so at around the I.R. fixed point. Finally, this problem is obviously present in effective theories where the cut-off 
reaches physical scales. 



B. Resumming the perturbation expansion 

The basic idea leading to the WH equation is the successive elimination of the degrees of freedom. This procedure 
produces the blocked action, the integrand of the functional integral of the blocked theory. This integrand must be 
well defined, ie each mode must either be left intact or be completely eliminated during the blocking. This is the 
positive side of the sharp cut-off in momentum space. The negative side is that it generates non-local interactions and 
spoils the gradient expansion. The simplest cure, to smear out the regulator and to use smooth cut-off procedure, is 
not a valid alternative. This is because the blocking based on smooth cut-off suppresses the modes partially and the 
integrand for the functional integral of the blocked theory is ill defined. Despite this general sounding argument one 
can proceed and generalize the successive elimination process for smooth cut-off in a rather surprising manner. 



1. Polchinski equation 

Let us start with the partition function 

Zk=J i5[$]e-^*«^~'*-*^^[*l, (38) 

where the interaction functional corresponds to the U.V. cut-off k. We split the propagator and the field variable 
into IR and U.V. components, 

Gk = Gk-Ak + Gk, ^ = <i) + 4>, (39) 

with the intention that the fields (p and (p should propagate with Gk-Ak and Gk, respectively. In other words, the 
kinetic energy contribution to the action is supposed to be of the form {<i>G^^^j^(j)+ |^G^^(^)/2 when written in terms 

of (j) and 4>. 

The introduction of the fields (j) and (j) separated by a smooth cut-off makes both (pp and (j)p non-vanishing and 
appears as a double counting of the degrees of freedom. In order to ensure the proper integration measure we introduce 
a dummy field, in such a manner that a suitable, momentum dependent linear combination of (j)p and (j)p reproduces 
and ^r. 



(40) 



and decouple 



l<t>Gkl^k<f> + = \^G-,'^ + \^G-^'^, (41) 



with a freely chosen dummy propagator Gd- Owing to the second equation in (39) Ai^x = Ai^2 = 1- The condition 
(41) gives 



Gk-Ak M = A2,,{p)\(G-,^ U 1 1 

G^V \l A2,2(n)A U2,l(D) A2,2(D) 



(42) 
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whose solution is ^2,1(1^) = \J GDGklGk-A.kGk, ^2,2(1^) = \J GoGk-Ak/GkGk- The transformation (40) is non- 
singular so long as Gk ^ Gk-Ak, the propagations of (f> and (j) arc distinguishable. 

Wo arc finally in the position to describe the blocking procedure which consists of the following steps: 

1. Redouble the number of degrees of freedom, <I> ($, $). 

2. Rotate the new, dummy degrees of freedom into the dynamics in such manner that the rotated fields split the 
original field variable and follow the prescribed propagation. 

3. Integrate out (p. 

In order to complete step 3. we insert a field independent constant into the partition function, 

Zk = J £)[$]£)[|.]e-A*^^'*-5k*G,-^<t-is^[*]. (43) 

The redistribution of the degrees of freedom by the inverse of the transformation (40) yields 

Zk = J D[<l)]D[4>]e-^'^^^-^>''^-^^'^^'^-i^'^^'''+^\ (44) 

up to a constant. This motivates the introduction of the blocked action Sk-Aki^] defined as 

^-Sk-Ak[<t>] = j D[4)\e-^^^k^4'-^Si[4>+4>]_ (45) 

The higher loop corrections were suppressed in the WH equation by restricting the functional space of the blocking, 
the volume of the loop integration in momentum space. But there is another way to suppress the radiative corrections, 
to decrease the propagator. Such a suppression method is better suited for smooth cut-off where the U.V. field can 
not be constrained into a restricted, 'small' functional space. This strategy is implemented by requiring that in the 
so far arbitrary split (39) Gk{p) = 0{/^k). The obvious choice is 



Gk{p) = i^kdkGkip). 



(46) 



By assuming the convergence of the perturbation expansion we expand the action in cj) and perform the integration 
over (j). Since the new small parameter is the propagator we retain graphs where the propagator appears in a minimal 
number. By assuming that the saddle point is at most O (Afc) we can expand the action in 0, 



^ + + ■^-itrlog(^G,-i + ^^J+0(Afe^) 

l^-9*=G,.^^-|tr(^a,Gfe^^jJ+0(Afe^) 



1 ;^TT ■ ^kGk 



e-iS.[^l+0(Afc2 



(47) 



The third equation was obtained by expanding the logarithmic function and gives the differential equation [9], [10] 



dkGk- 



The fourth line leads to a linear equation, 



fca.e-*^^[^] = nBe-i'i^'^\ S = 1 A . kdkGk ■ ^. 

2 00 0(p 



(48) 



(49) 



It is instructive to notice the similarities and differences with the WH equation. Eq. (45) looks as the starting 
point for the derivation of the WH equation and could have been obtained in a trivial manner for sharp cut-off. The 
crux of the argument leading to this relation is that this simple-looking equation is actually valid for smooth cut-off. 
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But in a surprising manner it is just the sharp cut-ofF hmit when the rest of the derivation of Polchinksi equation is 
invalid. This is because Gk as given in (46) is AfcO (AA:"^) = O (Afc°). 

Eq. (48) is a leading order equation in the perturbation rather than the loop expansion. It is formally similar to the 
WH equation, (10) when only the leading order is retained in the expansion (11) which corresponds to keeping graphs 
with a single vertex only in Fig. 1(b) and using the dashed lines to denote the propagation of (j). Notice that even if 
the higher order contributions of the perturbation expansion are suppressed in the limit Afc — > 0, the convergence of 
the perturbation expansion was assumed. 

The difference between WH equation and (48) is that the tree-level first term on the right hand side is always 
non- vanishing in Eq. (48), contrary to the case of the WH equation. It is easier to understand the reason of the 
non-trivial saddle point when the blocking (6) is considered in the presence of the source term. The source can 
be omitted and the cut-off independence of a single quantity, the partition function can be considered as the basic 
equation when the source has vanishing component in the functional subspace to integrate over in the blocking. In 
case of a smooth cut-off the background field has arbitrary momentum component and we can not really omit the 
source. One ought to follow the more cumbersome blocking which keeps the generator functional RG invariant. When 
classical physics (action) is modified by changing the cut-ofF, the field induced by the source changes, as well. This is 
the origin of the tree-level renormalization. But we made a short-cut and omitted the source. The price is that the 
tree-level renormalization can only be understood by letting the background field, (p vary freely. 

The understanding how the ^dependence may give rise the tree-level renormalization in (47) is an interesting 
demonstration of the well know theorem stating that the expansion in loops and in Ti are equivalent. In fact, the first 
equation in (47) introduces an effective theory for </> and the first two terms in the exponent of the second line give the 
efii'ective action corresponding to a given background field cj) in the classical approximation. In other words, the action 
with lowered cut-off which is supposed to reproduce the action with the original cut-off for arbitrary configuration 
(f) requires tree-level adjustment. The second term in the exponent has two factors of 'h~^ coming from the 'source' 
h~^5Sl[(l)\/5(() of the field and a factor of h from the </) propagator. The overall 1/h disappears in the differential 
equation (48) which displays Ti times the exponent. This is the way the graphs of Fig. 1(b) appear on the O 
tree-level in Eq. (48). Another reason is that they are graphs with no freely changing loop- momentum variable. One 
can understand in a similar manner that both the tree- and the loop-level graphs appear with a common factor of Ti 
inEq. (49). 



2. Gradient expansion 



We shall use here the results obtained in section II A 3 with the ansatz 

^1 



J X 



to derive the gradient expansion for Eq. (48). For the left hand side we find in the given order of the 
expansion 

(7(2) 



(50) 
gradient 



Z r Z -I f . • . t/(2) r 

- J V-pVpj^ + J V-pVpP ^P +VU+ — J V-pVp 

/ V-pVp- 



J V-pVp^ + VU + 



By means of the relation kdkG = —2Kp/k^ where = dK{z)/dz, z = p^ /k"^ the right hand side gives 



x§|i,,ot/l'>+>,-p(a(») + (Z-l)^ 



V—p—sVs 



C/(3) _ zW 



p - s {p + s) ■ s p- {p + s) 



K. 



Kp Kr 



Kr 



(51) 



14 




2 I " 



[/(4)+^(2)( P 



(52) 



This equation is consistent since the; terms O (77^) are comparable thanks to the smoothness of the cut-off. We compare 
the coefficients of the same ?7-dependent expressions and find 



and 



Jp \ 



(53) 



(54) 



We have two independent equations for two functions because the Eq. (54) follows from the first equation of (53). 
It is instructive to check the first equation in the local potential approximation Z — 1, 



U = 



|p=o 



[ kdkl 



(55) 



The first term describes the tree-level adjustment of the bare action in keeping the physics cut-off independent and 
some of the corresponding graphs are depicted in Fig. 1(a). The present form, characteristic of the gradient expansion, 
is non- vanishing when the change of the cut-off at p = A; is felt at the base point of the gradient expansion, p = 0. The 
second term is the leading order perturbative contribution to the WH equation (14) when the mass term is treated 
perturbatively, as expected and the beta functions. 



r kdk 

Jp p 



kdkKp 

,2 ' 



(56) 



for n > 2 correspond to the leading order first column in Eqs. (23) with 92 = 0. The independent mode approximation 
reproduces the leading order renormalization of the coupling constants and the numerical integration resums the 
perturbation expansion in the approximation where the effective vertices are replaced by their values at p = 0. The 
leading order perturbative contributions, retained in the RG equation are linear in the coupling constants and the 
classical, non-linear terms drop out in the linearized evolution around the Gaussian fixed point therefore the usual 
critical exponents are reproduced. 



C. Composite operator renormalization 

The renormalization of composite operators [11] seems to be a highly technical and subtle issue, dealing with the 
removal of the U.V. divergences from Green functions where composite operators are inserted. But it becomes more 
elementary and general [12] as soon as we are ready to give up the perturbative approach and can reformulate the 
procedure in a general, non-perturbative manner. In order to see this let us go back to the remark made at Eq. (6). It 
is not enough to impose the RG invariance of the partition function only since it yields a single equation for infinitely 
many coupling constants. Instead one should require the RG invariance of the generator functional for Green functions 
because any observable can be reconstructed from this functional. This turns out to be rather cumbersome due to the 
source dependence generated for the blocked action [13]. We took another direction by imposing the RG invariance 
of another functional, the blocked action, cf Eqs. (7) and (45). Instead of this strategy we shall now use the matching 
of observables computed at different values of the cut-off to rederive the RG equation. This is a non-perturbative, 
blocking inspired generalization of the multiplicative RG schemes mentioned in the Introduction in what the scale 
dependence of the observables is followed instead of those of non-physical bare coupling constants. 

We start with a toy model to demonstrate the natural relation between blocking and composite operator renormal- 
ization and the field theoretical application follows next by generalizing the source term in the generator functional. 
This term which usually involves the elementary field only is extended here for any local operator with inhomogeneous 
source and we consider it as part of the action which now contains inhomogeneous coupling constants. Our main 
point is that the beta functions corresponding to this extended action provide a bridge between the blocking and the 
composite operator renormalization. 
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1. Toy model 



Consider the two dimensional integral 

Z = j dxdj/e-^^^'^) (57) 

over X and y which play the role of low- and high-frequency variables, respectively, where the bare action is given by 
the expression 

11°° 
y) = -s^x^ + -syy"" + Y,9n{x + y)" (58) 

n=0 

in terms of bare coupling constants fif„. We use 7i = 1 in this section. The blocked action is defined as 



(60) 



The elementary, bare operators are {x + y)" with n = 0, 1, 2, . . . and their expectation values are 

Jdxdy{x + yYe-s^-^y'> _ J dx'^e-^^^^ 



The right hand side of this equation expresses the expectation value of the bare operator in terms of an operator of 
the blocked, 'thinner' system. The action S{x) can be expanded in the terms of the operators a;", 

S{x) = ls,x'+f2g'^x-^, (61) 



m=0 



giving rise the composite operators 



{^"} = ^ = J2 ^"'Sm,n (62) 



m=0 



where 



Sm.n = • (63) 

dgn 

We have chosen a basis for operators in the bare theory and searched for operators in the blocked theory. The 
opposite question, the construction of the composite operators [{x + y)"] of the bare theory which reproduce the 
expectation values of the blocked operators gives 

J dxdye-Si^'V) Jdxe-si^) ' ^ ' 

with 

oc 

[{x + y)"] =J2(^ + yr{S-')m,n. (65) 

m=0 

It is instructive to generalize the toy model for three variables, 

111°° 
S{x, y, z) = -s^x^ + -Syy"^ + -s^z^ + ^gn{x + y + zY, (66) 

with 
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e-S(x,, 



^'V) = J dze-^^=''y''\ e-^(^) = J dye-^^=''y\ (67) 
where the operator mixing looks hke 



dS{x,y) _ ^Idze-^(^'y'^) 
dg„ J dze-s^^'y'"^) ' 



These relations yield 



^ dgn ~ J dydze-Si-^V'^) ' ^^^> 



^ J dze-s(-'y'-) ' Jdye-S(-'y) ' ^ ^ 



indicating that the evolution of the operators comes from the elimination of the field variable in their definition. We 
can compute the expectation value of {x"} at any level, 



J dxdydze S{x,y,z) J dxdye ^('^■v) J dxe 



S(x) 



The lesson of this toy model is twofold. First, one sees that the formal definition of the renormalized operator 
at two different value of the cut-off differs because these operators are supposed to reproduce the same averages by 
means of different number of degrees of freedom. Second, it shows that the cut-off dependence of the renormalized 
operators can be obtained in a natural and simple manner by considering the bare coupling constants at one value of 
the cut-off as the functions of the coupling constants given at another value of the cut-off. 



2. Quantum Field Theory 

We generalize the operator mixing of the toy model for the scalar field theory given by the bare action 

SA[<t>] = / E Gn,.(A)0„(<^.) = J2 GniA)On{ct>.), (71) 

where On{<t>x) represents a complete set of local operators (functions of (px and its space-time derivatives) and Gn,x(A) 
denotes the coupling constant. We simplify the expressions by introducing a single index h for the pair (n, a;) labeling 
the basis elements for the local operators and = J^. The decomposition (p^ = 4'k,x + 4>k,x of the scalar field 
into a low- and high-frequency parts is carried out as in Eqs. (7) and (39). 
The Kadanoff- Wilson blocking 



g-Sfc[0fc 



j £>[^fe]e-^'^''^'=+'^'=l (72) 
for the action 

SM = Y.Gf^^^)On{<l>k) (73) 

h 

generates the RG flow 

kdkGn{k) = pn{k,G). (74) 

The blocked operators are defined by 

sn (A n ^^'^['^'^] /^WOn(0fc + ^fc)e-^^[^'°+^''l 
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in agreement with Eq. (68). They satisfy the equation 



J D[(j)k]e-S''^'l"'^ /£>[<^fc]D[^fe]e-'5A[<A.+<?.] ' ^ ' 

c.f. Eq. (70), showing that they are represented in the effective theory by the functional derivative of the blocked 
action with respect to the microscopical coupling constants. The form (73) of the action gives the operator mixing 

in a manner similar to Eq. (69). We introduced here the sensitivity matrix 

It is the measure of the sensitivity of an effective strength of interaction on the initial condition of the RG trajectory, 
imposed in the U.V. domain, c.f. Eq. (153). The composite operator {On{(j)k)}k introduced in Eq. (77) replaces the 
bare operator {On{(f))}A in the Green functions of the effective theory for the I.R. modes (f)k- 

The differential equation generating the operator mixing (77) is obtained in the following manner. The relation 

Ak 

Gn{k-M) = Gn{k)-—f3n{k,G) (79) 
is used to arrive at the sensitivity matrix. 



A A^ _ ^^'^(fc - ^fe) _ ^G^^k - Ak) 6Gi{k) 
A/c,AJ- ^^^^^^ -2^ ^^^^^^ ^^^^^^ 

■ _ Afc ,5fe(fc,G) - 
. ™^ k 5Gi{k) 



= E 



Slnik) (80) 



and to find 



kdkS^,n{k,K) = Y: A). (81) 

The scale dependence of the operator mixing matrix is governed by 

l ^fe(fc,G) 

= k^G-W ^ ^ 

Summary: Eq. (74) represents the link between Kadanoff- Wilson blocking and composite operator renormalization. 
On the one hand, Pn describes the evolution of the action in the traditional blocking scheme, Eq. (72). On the other 
hand, interpreting the coupling constants in the action as soiirces coupled to composite operators f3n determines 
the mixing of composite operators in Eq. (82). Therefore all conepts and results of the blocking, e.g. fixed point, 
universality, etc. has a counterpart in composite operator renormalization. 



3. Parallel transport 



The operator mixing discussed above has a nice geometrical interpretation, a parallel transport of operators along 
the RG trajectory [14] with the connection Tn^rh,{k)- The independence of the expectation value from the cut-off, 

/ D [<l>k' ] D [h' ] On ('/'fc' + ^fc' ) e- I-^^' 1 ^ / I3[0fe]g[0fc]O^((/)fc +0fc)e-^''[^''+'^''l 

defines the parallel transport of composite operators. The linearity of the mixing (77) assures that this parallel 
transport is indeed linear and can therefore be characterized by a covariant derivative 
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DkOk = {dk - T)Ok 



(84) 



where F^.n = {S~^ ■ 7 • S)n^m in such a manner that D^Ok = along the RG flow. 

The dynamical origin of the connection comes from the fact that there are two different sources the scale dependence 
of {Ok) comes from: from the explicit fc-dependence of the operator and from the implicit fc-dependence due to the 
cut-off in the path integration. The operator mixing is to balance them. In fact, the covariant derivative could have 
been introduced by the relation 

dk{Ok) = (DkOk), (85) 

requiring that the operator mixing generated by the connection is to make up the implicit fc-dependence of the 
expectation value coming from the cut-off. 

It is obvious that T is vanishing in the basis {On}k, 

dk{Y,Cn{k){On}k) = Y.dkCn{k){{On}k) = Cfi(fc){Oft}fc). (86) 

h n h 

The connection can in principle be found in any other basis by simple computation. 



4- Asymptotical scaling 



We linearize the beta-functions around a fixed point G^, 

/3n ~ ^ '^ImiGrh - G*^), (87) 

m 

and write the scaling coupling constant, the left eigenvectors of T 

E left p* _ ^ Q,-cl«*f! 

m 

as 

G^i = Y,ct%(Grn-G*^). (89) 

rh 

They display the scale dependence 

G^"~A:"". (90) 
Furthermore let us fix the overall scale of the local operators. This can be done by using the decomposition 

0(</'x) = ^&nO„(</.,) (91) 
n 

where 0„((^(.x)) is the product of the terms ■ ■ ■ dm(j)™{x) with coefficient 1. Thenorm||0|| = is introduced 

with the notation 

o = M- w 

for the operators of unit norm. We shall use the convention that the coupling constants Gft(A) always multiply 
operators of unit norm in the action. 
The scaling operators 

Or = E4t*0^ (93) 

rh 

are obtained by means of the right eigenvectors of F* , 
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Ep* right right /g.x 

m 

and they satisfying the conditions of completeness c'^'sht . ^left _ ^ ^^^^ orthonormality c'*^^* • c"^^^ = 1. The coupling 
constants of the action 

Sk = Y,Gn{k)Of{^ (95) 

n 

obviously follow (90). The operator 

= Y.^n'{OfU, (96) 

n 

written at scale k in this basis yields the parallel transport trajectory 

{0}u' = E bn{Jof}l}u' = E (^) (97) 

h h 

in the vicinity of the fixed point. 

The beta function introduced by Eqs. (72)- (74) agrees with the usual one. One finds that only relevant operators 
have non-vanishing parallel transport flow. Universality manifests itself at a given scaling regime in the suppression 
of the parallel transported irrelevant operators. 

The blocked action Sk[(l>k,x-,Gn{k),Gn{^)] possesses two interpretations: 

• The value of the coupling constant at the running cut-off, Gn{k), reflects the scale dependence of the physical 

parameters. 

• The dependence on the bare coupling constants, Gft(A), the initial condition of the RG flow provides us the 
generator functional for composite operators. 

Finally, the composite operator renormalization represents an alternative way to arrive at the RG equation. In fact, 
the beta functions arising from the blocking (72) are obtained in this scheme by the parallel transport, the matching 
(76) of the expectation values. 



5. Perturbative treatment 



Let us finally compare the matching of the obscrvablcs, described above with the traditional method of perturbative 
composite operator renormalization. The inversion of Eq. (77) gives 

[Oni^k + 4>k)]k = Y^Orh{^k + 4>k){S-\k,A))^,n. (98) 

rh 

It is not difficult to see that the operator [On{4>k + <Pk)]k of the bare theory which corresponds to the operator 
{On{(pk)}k of the eff'ective theory agrees with the result of the usual composite operator renormalization. In fact, 
the bare action (71) can be split into the sum of the renormalized part and the counterterms the framework of the 
renormalized perturbation expansion, giving Gub — GnR+GacT, Gn{A) = GnB and Gn{k) = GnR- The counterterms 
are introduced just to render certain Green functions with composite operator insertions cut-off independent [13]. The 
composite operator corresponding to the renormalized (ie blocked) operator On{(j>k) in this perturbative framework is 

^Orn{<Pk + 4'k)^-§P^ (99) 
m 

which agrees with (98). 

It is worthwhile recalling that the parameters and operators of a bare theory correspond to the cut-off scale in a 
natural manner. Therefore {On(0fe)}fc represents the bare operator [Ofi(0fe + <^fc)]fe in the effective description at the 

scale k and reproduces the observational scale dependence by construction in a manner similar to the scale dependence 
of the hadronic structure functions obtained in the framework of the composite operator renormalization convey the 
same information [15]. 
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D. Continuous evolution 



The RG equations obtained so far deal with the evolution of the bare action during the gradual lowering of the U.V. 
cut-off. There are two different reasons to look for alternative schemes where the evolution of the Green functions 
rather than the bare action is followed. 

One reason is usefulness. In the traditional RG strategy followed so far the renormalized trajectories give the 
bare coupling constants as functions of the rTinning cut-ofF. Though one can extract number of useful information 
from the trajectories they remain somehow qualitative because the bare parameters of the theory with the running 
cut-off are not physical quantities. Though this remark does not apply to the multiplicative RG scheme where the 
effective parameters are constructed by means of Green functions this scheme is seriously limited as mentioned in the 
Introduction. Returning to blocking, one can say at most that the difference between bare parameters and physical 
quantities arises due to the fluctuations in the path integral and this latter decreases with the number of degrees of 
freedom as we approach the IR end point of the trajectories. Therefore it is the blocked action with very low cut-off 
only which is supposed to be directly related to physical quantities (in the absence of IR instability). 

Another more formal point of view is to construct a strictly non-perturbative scheme. The RG equations (10), 
(48) represent the complete resummation of the perturbation expansion but they are not really non-perturbative 
equations. This situation is reminiscent of the Schwinger-Dyson equations. They were first obtained by resumming 
the perturbation c^xpansicm and only later by a genuine non-perturbative method, by the infinitesimal shift of the 
integral variable in the path integral formalism. The fact that the naive result derived by assuming the convergence of 
the perturbation expansion is correct is presumably related to the unique analytic continuation of the Green functions 
in the coupling constants. Can we find in a similar manner the truly non-perturbative RG equations? If possible, 
it will come by a different route: by relating two path integral averages instead of computing them by brute force. 
The quantities in question are one-particle irreducible (IPI) amplitudes and we shall give up to follow the evolution 
of the bare coupling constants, the parameters in the path integral. This will allows us to avoid any reference to the 
perturbation expansion. It remains to be seen if (10) and (48) can be derived in a similar manner or they remain 
valid for strongly coupled models. 

We begin this program by writing the generator functional for the connected Green functions in the form 



= j Li[,ji]e-^('5^['^1+'^^['^1--'-"^), (100) 

where the term Cfc[(/)] is introduced to suppress fluctuations. What we require is that (i) for k = oo the fluctuations 
be suppressed, Coo[0] = oo, (ii) the original generator functional be recovered for fc = 0, CQ[(j)] = and (iii) the 
fluctuations are suppressed only, ie C^; [(/)'] = for some configurations close to the vacuum expectation value, (f)' « (</)). 
The simplest choice for models without condensate is a quadratic functional, 

Ck[<iy] = \^-Ck-cl> (101) 

but in certain cases higher order terms in the field variables are necessary in the suppression. We distinguish two 
kinds of suppression. The length scale of the modes 'released' l{k) is well defined when 

The evolution generated by such suppressions shows the scale dependence in the spirit of the Kadanoff- Wilson blocking, 
the contribution of modes with a given length scale to the dynamics. Examples are, 

-P'M) [16] 
Cu{p)=\am' [17] (103) 



where /(p) = e~^'^'^^ 1^^^" and a,b,c> 0. Another kind of suppression for which Ml{p) shows no clear structure will 
be considered in section II E below. 

The evolution equation for W[j] is easy to obtain. 



S 



ei^^y]. (104) 
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This is already a closed functional differential equation but of little use. The reason is that it is difficult to truncate 
W[j] being a highly non-local functional. In order to arrive at a more local functional we shall make a Legendre 
transformation and introduce the effective action Tf. [0] as 



SW\j] 



with 



SW[j] 



—dkj + dkjcp = -dkWklj]- 



(105) 



(106) 



It is advantageous to separate the auxiliary suppression term Ck[(f)] off the effective action by the replacement r[(p] 
r[(/)] + Cfe[(/)] resulting in 



dkTk[^]=e-i'^''^i^dkCk e^"^"^^ - dkCk[(j>] 



(107) 



This relation takes a particularly simple form for the quadratic suppressions (101). Since the quadratic part of the 
effective action is the inverse connected propagator we find [16], [10], [19], [20], [21], 



dkTk[ 



tr [dkCk ■ m C 



n 



:tr 



9kCk 



Ck 



(108) 



Two remarks are in order in comparing the evolution equation obtained here with the pervious RG equations. 
First, the Kadanoff- Wilson blocking is constructed to preserve the generator functional W[j] for the thinner system 
and this is realized by eliminating modes at the U.V. side, by changing Sk[4>\- One is not aiming at keeping anything 
fixed in the continuous evolution schcnic. The IPI generator functional does depend on k which appears as the lowest 
momentum of modes considered in the theory whose effective action is Ffe [</>]. The second remark is that the right 
hand side is proportional to Ti. But the apparent absence of the mixing of tree and loop-levels is misleading since the 
propagator may be the sum of O {Ti^) tree-level and O (Ti) fluctuation contributions. The explicit tree-level term is 
missing because the modification of the bare action was carried out in the part O ((^^) and this correction to the free 
propagator was removed by the step T[(j)\ ^ r[^] -|- Ck[4>\, the substraction of the tree- level suppression term from the 
effective action. 

It is illuminating to compare the evolution and the Schwingcr-Dyson (SD) equations. Though the formal aspects 
appear similar the content of the equations differs. A similarity appearing immediately is that both set of equations 
determine the Green functions in a hierarchical manner. The derivation of the equations shows some similarity, as 
well. In fact, the crucial step in the derivation of the evolution equation is the first equation in (104). This step, 
the expression of the derivative of the functional integral formally, without its actual evaluation is the hallmark of 
the non-perturbative demonstration of the SD equations. Both the SD and the evolution equations are genuinely 
non-perturbativc because we actually do not evaluate the functional integral, instead wc relate to another one by 
bringing a derivative inside the path integral. In both schemes on compares two functional integrals which differ 
slightly, either in the infinitesimal shift of the integral variables or in the infinitesimal change in the action. Let us 
write the first equation in (104) for finite Afc, 



j D[4)] ^e-s(C'^[«^l+^-B[<^]-j-<?i) _g-^(Cfc_AfcM+SB[<^]-j-0)^ 
= \j D[ct>] (Afc5fcCfc[<^]e-*^'=[*l + O ((AA:)^)) e-*(^-l^ 



(109) 



The small parameter Afc is used to suppress the insertion of more field variables in the evolution equation, to cut off 
the higher order Green functions from the evolution, the strategy common with the non-perturbative proof of the SD 
equation. The higher order contributions in Afc bring in higher order Green functions and we suppress them by the 
smallness of the step in 'turning on' the fluctuations. This is how the RG idea, the replacement of the higher loops 
by running effective coupling constants and the dealing with small number of modes at each step, is realized in the 
evolution equation scheme. 

The obvious difference between the evohition and the SD equations is that the latter express the invariance of the 
functional integral in a manner similar to the RG equations but the former is simply an expression of the derivative 
of the functional integral with respect to a parameter. The similarity of the evolution equation and the RG strategy 
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leads to another difference. The subsequent elimination of modes and taking into account their dynamics by the 
introduction/modification of effective vertices generates a 'universal' evolution equations which does not depend on 
the theory in question. In fact, the evolution equation (108) and the functional RG equations Eqs. (6), (45) do not 
contain any reference to the model considered. The effective or the bare action of the model appear in the initial 
condition only. The SD equation is based on a careful application of the equation of motion within the expectation 
values and therefore contains the action of the model in an obvious manner. 

The evolution equation method seems to be better suited to numerical approximations then the SD equations. This 
is because the 'dressing', the summing up the interactions is achieved by integrating out differential equations rather 
then coupling the Green functions in the SD hierarchy, a numerical problem we can control easier. 

E. Blocking in the internal space 

We present now a generalization the RG method. The traditional RG strategy is aiming at the scale dependence of 
observables and consequently follows the cut-off dependence in the theories. The results are obtained by the successive 
modification of the cut-off and the accumulation of the resulting effects. One can generalize this method by replacing 
the cut-off by any continuous parameters of the theory on which the observables depend upon in a differentiable 
manner. Such a generalization replaces the RG equation by an evolution equation corresponding to the parameter or 
coupling constant in question. 

One gains and looses in the same time during such a generalization. One looses intuition of and insight into the 
dynamics since the trajectories generated by the evolution equation have no relation to scale dependence. To make 
things even more complicated the trajectories do not correspond anymore to a fixed physical content, instead they 
trace the dependence of the dynamics in the parameter considered. But we gain in flexibility. In fact, the parameter 
we select to evolve can be the Plank constant or a coupling constant and the integration of the resulting evolution 
resums the semiclassical or the perturbation expansion. We may apply this method in models with non-trivial saddle 
point structure or where the cut-off would break important symmetries, e.g. in gauge theories. 

A more technical aspect of this generalization touches the U.V. divergences in the theory. The traditional RG 
procedure is based on the tacit assumption that the U.V. divergences are properly regulated by the blocking, the 
moving cut-off. This seems to be a natural requirement if the RG flow is interpreted as scale dependence. In fact, an 
U.V. divergence left-over by the blocking would indicate the importance of the processes at the regulator and would 
lead to the appearance of a second scale. In the generalization of the RG procedure we introduce below there may 
not be a well deflned scale related to the evolution and we have to regulate the U.V. divergences. 

There are two different kind of physical spaces in Field Theory. The space-time or the momentum-energy space 
where the 'events' arc taking place will be called external space. The 'events' arc characterized by further quantities, 
the field amplitudes. The space where the field amplitudes belong will be called internal space. Therefore the field 
configuration realizes a map cp '■ external space — > internal space. The traditional Kadanoff- Wilson blocking 
orders the degrees of freedom to be eliminated according to their scale in the external space. We may realize a similar 
scheme by ordering the modes according to their scale in the internal space, their amplitude [22]. We present two 
different implementations of this idea, one for Wilsonian and another for the effective action. 

1. Wilsonian action 

We shall use the flexibility of the method of deriving the Polchinksi equation in section II B to separate off and 

eliminate the fluctuation modes with the largest amplitude in the path integral. Since the fluctuation amplitude is 
controlled by the mass this leads to the description of the dynamics by highly massive modes, small fluctuations at 
the final point of the evolution. 

Let us consider the model (38) where the role of the parameter k is played by the mass M and = p'^ + M^. 
We shall follow the evolution from = to >> where A is the U.V. cut-off. The decomposition 

Gm^ = Gm^+am^ + Gm2, ^ = (p + (j), (110) 

will be used with 

1 1 ~ fM^ — O)"^ ~ 1 1 

(M^ + AM'-n)cp+ -0L__^^ = (m2 - □) $ + (Ill) 

We replace the degrees of freedom <& by the field with infinitesimal fluctuations. After eliminating (j) the remaining 
field has smaller fluctuations than We use the ansatz 
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-{ZM[(j)x) - '^)dfjt4>xd^i4>x + UM{(t>.,:[ 



(112) 



where the U. V. regulator with smooth cut-off is not shown exphcitly. We follow the steps outlined in section II B 2 
which gives rise the evolution equation with left hand side 



dM-{lJ V-pVp [{Z - ly + + VU^ . 



The right hand side reads as 
1 



J 5,,oU^') + Vp [{Z - ly + C/(^)] +IJ^ VrVp-r [U^'^ - Z^') {vp-T^+ 



(113) 



(p2 + 



xUp,o;7W+r?-p [(Z - l)p2 + [/(2) 



+ 2 / '^-P-s^s 



C/(3) + ^(1) ^s^^p-S+p" 



2 (p2 + M2)2 



V 



\Z - l)p2 + [7(2 



2 / 'n-rVr 



[/(4)+^(2)(p2^^2) 



1 



{Z - l)p2 + [/( 



2) 



^ / '7-r?7r 



[/(4)+^(2)(p2^^2 



2 7p (p2 + M2)2 

The system of evolution equation projected onto the different jy-dependent terms is 

(Z - + [/(2) 



2M^0..U = -U^^^^ + nM^l (p2 + M2)2 

7"^(2)2 r 1 

^ 2f/Wz(i)+?iM4z(2) / 

Jj, (p2 + M2)2 



2M'^dM^Z = -AU^'^\Z - 1) + 4- 



M2 



ip (p2 + M2) 



(114) 



(115) 



The O (?i") terms on the right hand side are to keep the tree- level observables M-invariant within the gradient 
expansion ansatz. As far as the loop corrections are concerned, the beta functions 



/3„ = MdM9n = nM^U^-+-) ^ p^M^' 



(116) 



obtained from the first equation with Z = 1 agree with the leading order contribution to the one-loop renormalized 
potential, 



(117) 



except their sign. This is because we intend to keep the partition function unchanged as opposed to Eq. (117) where 
the mass is evolving together with the dynamics. 



2. Effective action 

Let us consider for the sake of simplicity again the scalar model given by the bare action 
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Sb = 



(118) 



and the suppression with no structure in the external space [22] , 

Ck{p) = -j- 



(119) 



which acts as a 'smooth cut-off' for the ampUtude of the fluctuations. The corresponding evolution equation is 



Ml 



(120) 



The effective action of the theory is obtained by integrating this equation from the initial condition Tmo = '^b [4>] 
imposed at M"^ = » » m| to M ^ 0. 



Let us consider how this scheme looks like for the ansatz 



rM['A] 



(121) 



One finds 



= -[ 
2 J J, 



^ = 2 



(Zp2 + M2 + C/(2))^ (Zp2 + M2 + m2))^ 

(Z(V + £/(3))^ _ 8p2_,,, (z(^y + t/(3)) 

(Z392 + M2 + C/(2))* d {Zp^ + M'^ + U(^))^ 

d (Zp2 + M2 + C/(2))^ 

The asymptotic form of the first equation for » C/^^^ in the local potential approximation, 

M2 



MOmUm « -U 



I 



(p2 + M2)2 



(122) 



(123) 



up to field independent constant and it agrees with asymptotic form of the first equation in (115) except the sign. 
The integrand in (123) corresponds to transformation rule (4) of the propagator under the infinitesimal change of the 
mass for the first, leading order graphs in Figs. 3. The numerical integration resums the higher orders in the loop 
expansion. In short, we see the Callan-Symanzik scheme at work in the functional formalism. 

The suppression (119) freezes out modes with momentum p if p'^ < and generates a characteristic scale Pcr 
where {Zp'^ + M ^ + C/^^) ) / (^Zp"^ + J/^^) ) deviates from 1. So long this scale is far from the intrinsic scales of the model, 
>> Z7^2) in U.V. scaling regime, such an internal space blocking generates the external scale Pcr ~ M/y/Z. In 
other words, the universal part of the beta functions obtained in the U.V. regime by means of the Callan-Symanzik 
scheme should agree with the same part of the beta functions coming from other schemes. This agreement has been 
observed in the framework of the scalar model [23] . We shall check this quickly in the local potential approximation. 
By comparing Eq. (123) with the asymptotic form of the WH equation (14), 



kdkUk = - ^ 



(2) 



2(27r)« 



we find 



A/k 



-^d,y 



(124) 



(125) 



Below the upper critical dimension, d < 4, the right hand side is finite and the external (WH) and the internal 
(CS) scales are proportional. At the critical dimension the M dependence of the right hand side is through the 
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proportionality factor ln(l + A^/M^). The two scales become approximately proportional only for much higher 
value of the cut-off A, since the RG flow spends more 'time' close to the non-universal short distance regime due to 
the tree- level marginality of g^. There is no agreement beyond d = 4 where the cut-off scale is always important. 
The agreement between the two schemes is violated by higher order terms in the perturbation expansion since they 
represent the insertion of irrelevant effective vertices and by leaving the asymptotical U.V. regimes. This explains 
that similar argument does not hold when Z displays important M-dependcncc. 

It is worthwhile noting that the non-trivial wave function renormalization Z and the corresponding anomalous 
dimension r] of the field variable can be thought as the reflection of the mismatch between the scale dependence in 
the external and internal spaces around the U.V. fixed point. In fact, the phenomenological form 

{cl>^cj>y) « c\x - t/l^-''-'' (126) 

cf Eq. (128) below connects the fundamental dimensional objects of the internal and external spaces. 

The flexibility of choosing the suppression functional Ck [(j)] may be important for certain models. By choosing 

cm = jS[<l>] (127) 

where A is the U.V. cut-off the evolution in k resums the loop expansion since the effective Planck-constant. h^^{k) — 
h^^ +k/A evolves from to h. This scheme is advantageous for models with inhomogeneous saddle point eg solitons 
or instantons because their space-time structure is 'RG invariant', being independent of the gradual control of the 
amplitude of the fluctuations. Another advantage offered by the flexibility in choosing the suppression functional is 
the possibility of preserving symmetries, the point considered in section III E below. 

III. APPLICATIONS 

We shall briefly review a few applications of the functional evolution equations. An incomplete list of the develop- 
ments not followed due to limitation in time is the following. The exciting competition for the 'best' critical exponents 
have led to several works using this method [24]. [53], [17], [50], [51], [25], [57], [26]. Such kind of application opens 
up the issue of understanding the impact of truncation on the blocking [27], [28], [29] [25], [30], [31] and looking for 
optimization [32], [33], [18]. The phase structure and the nature of phase transitions are natural subjects [21], [34], 
[35] . The incorporation of fermions is essential to arrive at realistic models in High Energy [36] and Condensed Matter 
Physics [37]. The computation of the quenched average of Green functions, one of the essential obstacle of progress in 
Condensed Matter Physics can be approached in a new fashion by applying the internal space renormalization group 
method [38] . A promising application is in general relativity [39] . Finally, the infamous problem of bound states can 
be reconsidered in this framework [40], [41]. 

Much more to be found in review articles [42] , conference proceedings [43] and PhD thesis [44] . 

A. Fixed points 

1. Resettling 

'The original form of the RG procedure consists of two steps, the blocking, followed by a rescaling. The latter can 
be omitted if the RG strategy is used to solve models only. But it becomes important when we try to understand the 
scale dependence of the theories. Our physical intuition is based on classical physics and particles. After losing much 
of the substance during the quantization procedure there is a chance to recover some clarity by the introduction of 
quasiparticles, localized excitations with weak residual interaction. The rescaling step of the RG scheme is to check a 
given quasiparticle assumption by removing the corresponding scale dependence. The remaining scale dependence is 
a measure of the quality of the quasiparticle picture. The deviation from this non- interactive system is parametrized 
in terms of anomalous dimensions. 

The particle content of a theory is usually fixed by the quadratic part of its action. Therefore the rescaling is defined 
in such a manner that this part of the action stays invariant during the blocking. For an 0{d) invariant Euclidean 
or relativistically invariant real-time system the resulting rescaling reflects the classical dimension of the dynamical 
variables and coupling constants. For non-rclativistic systems the fixed point scaling is an artificial device only to 
identify the non-interacting part of the action. When smooth cut-off is used its details may influence of the scaling 
properties of the quadratic action, as well, and may induce deviations from classical dimensions.. 
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Consider the scale transformation 

p^p' = {l + e)p, x' = {l-e)x, + e)x') ^ {1 - ed^)(l){x) (128) 

in an 0{d) invariant scalar model where = (d — 2 + r])/2. The parameter r] reflects the deviation of from classical 
dimensional analysis. The efl'ect of this rescaling on the action 

S[4>\=^ Upi,-,Pn^Pl + -+Pn,0(i>Pl---(t>Pn^ (129) 

„ Jpi,-,Pn 

can be found in the following manner [45] : 

• The momentum integral measure changes as d'^p — > (1 — de)d'^p', giving 5 — > [1 — ^d (l)p-^]S . 

• The coupling constants change as — > W(i_£)p'^^.,, This can be written as the transformation 
S ^ [I — e /p 4'pP • d'p-^-]S of the action where the prime on the gradient dp indicates that the derivative acts 
on the coupling constants only and not on the Dirac-deltas. 

• The Dirac-deltas change as 5p^^ \-p„,o ~* ^{i-e)p[+—+{i-e)p'„,Oj amounting to the transformation S* — » [1 + de]S 

of the action. 

• Finally, the field transforms as (j)p ^ [1 — e{d^ — d)](j)p, inducing 5* ^ [1 — e{d^ — d)] (pp-^S. 
Adding up these contributions we find the rescaling generator 

Jp 

The RG equation for the rescaled action can be obtained by adding Q acting on the (effective) action to the right 
hand side, eg the evolution (49) turns out to be 

fcafee-^^^I-^l = (g + hB)e-i^'^^'l'\ (131) 

The rescaling of the field may be viewed as a transformation of the action without changing the physics at the fixed 
point. The terms generated on the fixed point action by the rescaling, gS*[4>\ are called redundant at the fixed point 
in question. 



S S 

^ d<pp d<pp 



(130) 



2. Reparametrizing 



We considered linear rescaling of the field variable but one can easily generalize the rescaling to non-linear 
reparametrization. The infinitesimal change (j)x ^ (l>'x = (l>x + ^^[^^ A performed inside the path integral gives 



D[ 



5'i'[(j);x] 



-S[0]-6/^*[0;x]^ 



The reparametrization invariance of the integral assures that the modification 



S[cj^] ^ S[<P] + eg^S[cP] 



of the action where 



g<s/ 



(132) 



(133) 



(134) 



is an invariance of the partition function. It has furthermore been noted [46] that some of the infinitesimal blocking 
relations for the (effective) action can be written in the form (133). 

Though it is certainly very interesting to find a common structure for the functional RG equations the connection 
with reparametrization invariance is not clear. The point is that non-linear rcparamctrizations can not be carried out 
inside of the path integral as in ordinary integrals. The problem arises from the fact that the typical configurations 
are rather singular in the path integration. Calculus known from classical analysis is replaced by Ito-calculus in 
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Quantum Mechanics due to the nowhere-differentiable nature of the quantum trajectories [47]. The quantum field 
configurations are even more singular. Consider a free massless field (p in d dimensions and its partition function 

Z = l[[ #xe-^ = JJ / dct>.,-i E.(A.^.)= (135) 

where lS.^(f>x = (t>x ^ (f'x-n and the dimension of the field was removed in the second equation, (j) — a'^/^~^(j). The 
typical configurations have A(j) — O (oP), ie the discontinuity of the original field variable is /S.(f) — O (a^^''/^). 

This simple scaling, the basis of the usual U.V. divergences in quantum field theory, gives non-differentiability 
in Quantum Mechanics, for d = 1, finite discontinuity in d = 2 and diverging discontinuities for d > 2. The non- 
differentiability can be represented by effective potential vertices in Quantum Mechanics, [48], a refiection of the 
unusual quantization rules in polar coordinates, sensitivity for operator ordering and quantum anomalies [49] . In two 
dimensions the discontinuity is finite and the continuous structure can either be preserved or destroyed by quantum 
fluctuations, cf. section HID. In higher dimensions the singularities remain always present. It is this singular nature 
of the trajectories which requires that one goes to unusually high order in the small parameter e characterizing the 
infinitesimal change of variables in the path integral. This problem, composite operator renormalization, renders the 
non-linear change of variables a poorly controlled subject in quantum field theories. 



3. Local potential approximation 

We shall search for the fixed points of the Polchinski equation in the local potential approximation, [50] [51]. We 
introduce dimensionless quantities (j) = k^~'^/'^^ and Uk{4>) = k'^Uki^)- The RG equation with rescaling for u is 

u = -du+ !1±^Z1^W _ „(i)2j^^ + f^j^'u<-^^ (136) 

where K' = k~'^ J^K'^ and we set = in the approximation Z = 1. Notice the explicit dependence on the cut- 
off function K. This refiects the fact that the RG flow depends on the choice blocking transformation. Only the 
qualitative, topological features of the renormalized trajectory and the critical exponents around a fixed point are 
invariant under non-singular redefinitions of the action. 

Gaussian fixed point: The fixed point equation, u* = has two trivial solutions, 

u*i(P)=0, (137) 

and 

where Kq, K' < 0. For any other solution u = O (^'^) as </> ^ oo. 

In order to find the scaling operators we introduce a perturbation around the fixed point by writing u = u*+ek~^v{(f)) 
where e is infinitesimal and solve the linearized eigenvalue equation 



(139) 



Having a second order differential equation one can construct a one-parameter family of solution after having imposed 
say v^^\0) = 1 in a theory with the symmetry (f> — > —(j). But the polynomial solutions 

n 

^'n(<^) =^i'M'A'^ (140) 

which are parametrized by a discrete index n correspond to discrete spectrum [52] . The critical exponents are identified 
by comparing the terms O ((t>^") in Eq. (139), A„ = d + {2 — d)n, and A„ = d — (rf + 2)n, for the fixed points (137) 
and (138), respectively [53]. The leading order critical exponents are given entirely by the tree-level contributions. 
The dimensionful coupling constants are cut-off independent in this case and we have An = [g2n] for the fixed point 
(137) where the dimension of the coupling constant [g] is given by Eq. (24). The dominant, largest exponent At and 
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the corresponding coupling constant play distinguished role in the scaling laws. The dimcnsionlcss coupling constant 
can be identified with the 'reduced temperature' t. According to the definition ^ w t^'^ of the critical exponent i/ we 
have the mean-field exponent, v — 1/Xt = 1/2. The mass is relevant at the point (137) as expected, as shown by the 
beta function (32 = —'^K^oi " ^52, according to Eq. (55). The scaling potentials with exponential growth for large (j) 
[54] correspond to the continuous spectrum and their physical interpretation is unclear. 
All exponents are negative at (138), this is an IR fixed point. 

Wilson-Fischer fixed point: There are non-Gaussian fixed point solutions for Eq. (136) when 2 < d < 4. One can 
construct a one-parameter family of solutions but the fixed point potentials which remain non-singular for arbitrary (f) 
correspond to a discrete set [55] , [52] : the Wilson-Fischer fixed point for 3 < d < 4 and as many fixed points as relevant 
when 2 < d < 3. The perturbation around these fixed points is a one-parameter family of scaling potentials with 
continuous spectrum of critical exponents. The restriction for solution which are finite and non-singular everywhere 
produces a discrete spectrum [56] in good agreement with other methods of determining the critical exponents [53], 
[57]. 

The truncation of the fixed point solution at a finite order of </) introduces error and spurious solutions [27], [28] 
which can partially be eliminated by expanding the potential along the cut-off dependent minimum [29]. 



4- Anomalous dimension 
The RG equation for the wave function renormalization constant z{(j)) = Z{'^) is 

z = -4K',u^^\z - 1) - 2i^^u(i)z(i) - 2&^u^^^^ 

fv 

+nz^^'>K' - r,{z - 1) + (141) 

The lesson of the case Z = \ \s that the requirement of the existence and finitcncss of u*{(f)) introduces discrete 
number of fixed point solutions. Let us try to follow the same strategy again. The dominant terms of such fixed point 
solutions w = i = of Eqs. (136) and (141) are 

u*{^)^'^^(t>'^+A(t)'^^, z*{(l>)^B. (142) 

Together with the conditions ^^^^(0) = ^(^^(O) = imposed for the ^ —cj) symmetrical models the solution are well 
determined in terms of A and B. The problem is that such kind of argument does not fix the value of 77. 

There is another condition to be fulfilled by the fixed points, the critical exponents should be invariant under 
rescaling. This is sufficient to determine -q. Unfortunately the rescaling invariance of the fixed point is lost unless 
sharp or specially chosen polynomial smooth cut-off is used [56], [58]. Furthermore the truncation of the gradient 
expansion contributes to the violation of this invariance, too. An approximation to find the 'best' solution when 
the rescaling invariance is not respected is the following [51]: Introduce a further condition which violates rescaling 
invariance, say fix the value of z{Q). This allows the determination of rj which would be unique if rescaling could 
be used to relax our last condition. We are as close as possible to the invariant situation within our parametrized 
problem when the dependence of r] on ^;(0) is the slowest. Therefore the condition dr]/dz{0) = selects the 'best' 
estimate of rj. 

The numerical error due to the truncation of the gradient expansion has been the subject of extensive studies in the 
framework of the scalar model, c.f. Ref. [42]. The general trend is that the truncation of the gradient expansion at the 
local potential approximation (zeroth-order) or at the wave function renormalization constant level (first order) yields 
approximately 10% or 3% difference in the critical exponents compared with Monte-Carlo simulations, seven-loop 
computations in fixed dimensions or the fifth-order results of the expilon-espansion. 



B. Global RG 



The usual application of the RG method can be called local because as in the determination of the critical exponent 
it is performed around a point in the space of coupling constants. Models with more than one scale may visit several 
scaling regimes as the observational scale changes between the U.V. and the IR fixed points. The determination of 
the 'important' coupling constants of such models which parametrize the physics goes beyond the onc-by-one, local 
analysis of the scaling regimes. It requires the careful study of crossovers, the overlap between the scaling laws of 
different fixed points, a problem considered in this section. 
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1. RG in Statistical and High Energy Physics 



It is important to realize the similarity the way RG is used in Statistical and High Energy Physics despite the 
superficial differences. The most obvious difference is that while the running coupling constants are introduced in 
Particle Physics by means of Green functions or scattering amplitudes at a certain scale, the parameters of Solid 
State Physics models are defined at the cut-off, the latter being a finite scale parameter, eg lattice spacing. Since the 
bare coupling constants characterize the strength of the physical processes at the scale of the cut-off, the two ways of 
defining the scale dependence are qualitatively similar and reproduce the same universal scaling laws. 

The U.V. fixed point where the correlation length diverges in units of the lattice spacing in Statistical Mechanics 
corresponds to renormalized theory in High Energy Physics where the U.V. cut-off is removed. 

There are two classification schemes of operators, one comes from Statistical and the other from High Energy Physics: 
if the coupling constant of an operator increases, stays constant or decreases along the renormalized trajectory towards 
the IR the operator is called relevant, marginal or irrelevant in Statistical Physics. The coupling constants which can 
be renormalized as the U.V. cut-off is removed in such a manner that observables converge are called renormalizable 
in High Energy Physics. The important point is these classification schemes are equivalent, in particular the set of 
irrelevant operators at an U.V. fixed point agrees with the set of non-renormalized ones. 

One can easily give a simple quantitative argument in the leading order of the perturbation expansion. First let us 
establish the power counting argument about renormalizability. Suppose for the sake of simplicity that there is only 
one coupling constant, g, playing the role of small parameter and a physical quantity is obtained as 

(0)=^5"/n (143) 

n 

where In is sum of loop integrals. This series gives rise the expression 

[/„] = [O] - n{g\ (144) 

for the (mass) dimension of the loop integral. Let us recall that the degree of the overall (U.V.) divergence of a loop 
integral is given by its dimension. We distinguish the following cases: 

• [(?]< : The higher powers of g decrease the dimension in (143) which is compensated for by increasing the 
degree of the overall divergence of the loop integrals, cf Eq. (144). Graphs with arbitrary high degree of 
divergence appear in the perturbation expansion and g is called non-renormalizable. 

• 5 [.9] > : The higher order loop integrals are less divergent, there are finite number of divergent graphs in the 
perturbation series of any observable. The coupling constant is super-renormalizable. 

• [51] = : The maximal degree of (U.V.) divergence is finite for any observable but there are infinitely many U.V. 
divergent graphs. 5 is a renormalizable coupling constants. 

Since the IR degree of divergence of massless loop integrals is just — the perturbation expansion of a super- 
renormalizable or non-renormalizable massless model is IR unstable (divergent) or stable (finite), respectively. The 
compromise between the U.V. and the IR behaviors is attained by dimensionless, renormalizable coupling con- 
stants. The IR stability can be realized by a partial resummation of the perturbation expansion in massless super- 
renormalizable models [59]. 

It is rather cumbersome but possible to prove by induction that the definition of renormalizability outlined above 
according to the overall divergence of the loop integrals remains valid in every order of the perturbation expansion. 

In order to separate off the trivial scale dependence one usually removes the classical dimension of the coupling 
constants by means of the cut-off in Statistical Physics. On the tree- level, in the absence of fluctuations classical di- 
mensional analysis applies giving the relation g = k^^^g between the dimensional and dimensionless coupling constants, 
g and g, respectively. Assuming that there is no evolution on the tree-level we have the scaling law 

/ u\-[9]+0{h) 

m = [j) 5(A) (145) 

for the dimensionless coiipling constants, showing that the non-renormalizability of a coupling constant, [g] < 
is indeed equivalent to irrelevance, the decrease of the coupling constant towards the IR direction. Higher loop 
contributions do not change this conclusion so long as the loop expansion is convergent and the anomalous dimension, 
the contribution O (h) in the exponent can not overturn the sign of the classical dimension [g] . In case of a marginal 
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coupling constant, [g] = on has to carry on with the loop expansion until the first non- vanishing contribution to the 

anomalous dimension. 

One can construct a much simpler and powerful argument for the equivalence of the irrelevant and the non- 
renormalizable set of coupling constants in the following manner. Consider the U.V. fixed point P and the region 
around it where the blocking relations can be diagonalized, as depicted in Fig. 5. The solid line shows a renormalized 
trajectory of a model which contains relevant operator only at the c\it-ofF. The Lagrangian of another model whose 
trajectory is shown by the dotted line contains an irrelevant operator, as well. The difference between the two models 
becomes small as we move towards the IR direction and the physics around the end of the U.V. scaling regime is more 
independent on the initial value of the irrelevant coupling constants^ longer the U.V. regime is. This is what is called 
universality of the long range, low energy phenomena. By looking "backwards" and increasing the cut-off energy, 
as done in Particle Physics, the non-vanishing irrelevant coupling constants explode and the trajectory is deflected 
from the fixed point. As a result we cannot maneuver ourselves into the U.V. fixed point in the presence of irrelevant 
operators in the theory. Since the infinite value of the cut-off corresponds to the renormalized theory, represented by 
the fixed point, the irrelevant operators are non-renormalizable. Despite its simplicity, this argument is valid without 
invoking any small parameter to use in the construction of the perturbation expansion. 




FIG. 5. An U.V. fixed point and its vicinity. The x and the y axis correspond a relevant and an irrelevant operator, 
respectively. The circle denotes the vicinity of the fixed point where the blocking relation is linearizable. 

Renormalizable models were thought to be essential in High Energy Physics. These models in which a sufficiently 

large class of observables can be made convergent are distinguished by their simplicity, their interactions exist without 
any cut-off. Examples are Yang-Mills models in four dimensions and the sine-Gordon, Thirring and the X-Y models in 
two dimensions. The non-renormalizable models are those whose existence requires a finite value of the U.V. cut-off, 
like QED, the 0^ model in four dimensions, the Standard Model and models in Solid State Physics. It is believed 
that asymptotically free models are renormalizable only. 

The traditional reason to discard non-renormalizable models was their weak predictive power due to the infinitely 
many renormalization conditions they require. Since we can never be sure what kind of Lagrangian to extrapolate 
up to high energies we need another point of view. According to the universality scenario the non-renormalizable 
models can be excluded because they cannot produce anything different from renormalizable theories. The cut-off 
independent physics of the latter is parametrized by means of renormalized coupling constants. But the subtle point 
to study below is that this reasoning assumes the presence of a single scaling regime with non-trivial scaling laws in 
the theory which is a rather unrealistic feature [7]. 

Thus renormalizability is the requirement of a simple extrapolation to the U.V. regime without encountering " new 
physics" . The evolution of High Energy Physics shows that this is a rather unrealistic assumption, any model with 
such a feature can be an approximation at best, to study a given interaction and to sacriiy the rest for the sake of 
simplicity. The goals are less ambitious in Statistical Physics and apart of the case of second order phase transitions 
the renormalizability of models is required seldom, for convenience only, not to carry the regulator through the 
computation. All models in Solid State Physics are effective ones given with a physically motivated cut-off. 

Let us start the discussion of the possible effect of the co-existence of several scaling regimes along the KG fiow with 
the simplest case, a model with a gap in the excitation spectrum above the ground state, ie with finite correlation 
length, ^ < oo. Suppose that the RG fiow of the model starts in the vicinity of an U.V. fixed point and reaches an 
IR fixed point region as shown in Fig. 6. Universality, the parametrizability of the physics beyond the U.V. scaling 
regime by the relevant coupling constants of the U.V. scaling laws, tacitly assumes the absence of any new relevant 



'^More precisely all the irrelevant coupling constants can modify is an overall scale. 
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coupling constants as we move along the renormalized trajectory towards the IR. as shown in Fig. 5. This assumption 
is in fact correct for massive models whose IR scaling laws are trivial, meaning that that only the Gaussian mass 
term is relevant. To see why let us consider the RG flow in the IR side of the crossover, when the U.V. cut-off, say 
the lattice spacing a >> ^ and make a blocking step, a a' which generates the change gn{a) gnW) in the 
coupling constants. The change A(7„(a) = 5n(o') — gn{a) is due to fluctuation modes whose characteristic length is 
a < £ < a'. Owing to the inequality ^ << a < £ these fluctuations arc suppressed by e~^^^ and the flow slows down, 
gn{a) ~ gn{a')- This implies the absence of run-away trajectories, the absence of relevant non-Gaussian operator. 
But more realistic models with massless particles or with dynamical or spontaneous symmetry breaking occurring at 
finite or infinite length scales, respectively, or with condensate in the ground state the IR scaling may generate new 
relevant operators and universality, as stated in the introduction, referring to a single scaling regime is not a useful 
concept any more. 



uv 




IR 



FIG. 6. The RG flow of a model with massive particles. The renormalized trajectory connects the U.V. and IR scaling 
regimes, indicated by circles which are separated by a crossover at the intrinsic scale of the model, p « 1/^. 



2. Scalar model 



In order to understand better the generic case shown in Fig. 6 we return to its simplest realization, the scalar 
model in the local potential approximation, Eq. (14), [7]. One does not expect anything unusual here since the model 
possesses a gap and the IR scaling regime is trivial in cither the symmetrical or symmetry broken phase. But wc shall 
see that the model has soft large amplitude fluctuations and the true vacuum is just there where this non-perturbative 
and the well known perturbative domains join. The second derivative of the local potential is discontinuous there, 
depends on the direction from which the vacuum is approached. We shall arrive at some understanding of the 
soft modes in two step, by separating the loop contributions from the tree-level structure. We start with the more 
traditional loop contributions. 

To facilitate the dealing with soft modes we introduce an external constraint which controls the expectation value 
Jxi4'x)j = V^. The usual realization of this constraint is the introduction of an external source as in Eq. (5). 
We shall consider the model in the symmetry broken phase where J^{4>x)q = ^^vac 7^ and choose homogeneous 
external source jx = J which can destabilize the naive symmetry broken vacuum whenever J^vac < and can induce 
1^1 < $vac- Note that such a control of the field expectation value by an external source preserves the locality of the 
action. 

The beta-functions, as given by Eq. (23), depend strongly on the propagator G{k'^) = l/{k^ + g2{k)), represented 
by the the internal lines in Figs. 3. The qualitative features of the function G~^{k^) are shown in Fig. 7 for $ = 0. 
According to the loop expansion 



(146) 



the bending of the lines in the figure is due to radiative corrections, as far as the perturbative region, » 1 is 
concerned. The dashed line corresponds to the massless Coleman- Weinberg case [23]. The flow above or below this 
separatrix is in the symmetrical or symmetry broken phase, respectively. This qualitative picture is valid so long 

1*1 < *vac- 

It is important to recall that the argument of the logarithmic function in the WH equation in Eq. (14), G~^{k'^), 
is the curvature of the action at the cut-off and serves as a measure of the restoring force driving fluctuations back to 
their trivial equilibrium position. Thus perturbative treatment is reliable and the loop contributions are calculable so 
long G-\e) = G-i(fc2)/A;2 » 
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In order to have a qualitative ides if what is happening let us consider the model where Ua{^) = m^0^/2 + g(j>'^/A\ 
at the cut-off, fc = A. By simply ignoring the loop corrections we have g2{k) = m^, 9A{k) = g and gn{k) = for 
n > 4 on the tree-level, ie and g can be identified with renormalized parameters and < 0. For an isolated 
system (j = 0) where 1$! = $vac = —6m^/g and the use of the perturbation expansion appears to be justified for 
G~^(fc^) = 1 — w? /k"^ > gik"^^^. Spontaneous symmetry breaking, or the appearance of a condensate in the ground 
state in general, is characterized by instability of fluctuations around the trivial, vanishing saddle point. When an 
external source is coupled to 0x to dial |$| < 3>vac in the symmetry broken phase then non-perturbative effects set in 
as k decreases. We can see this easily when the field expectation value is squeezed in the concave part of the potential, 
1^1 < $infl = —2m?/g = $vac/3. In fact, approaches zero and the amplitude of the fluctuations explodes at 
k = sj—rn?- — g(^'^/2. Below this non-perturbative regime the plane wave modes become unstable since G~^{k'^) < 
and a co herent state is formed, reflected in the appearance of a non-homogeneous saddle point with characteristic 
scale l/-^-m^ — g$^/2 

The mechanism responsible of spreading the instability over the whole region |$| < $vac comes from the tree-level 
structure and will be discussed later, in section III C. 




FIG. 7. The inverse propagator G^^(p^) of the scalar model as the function of the momentum square for "1> = 0. The 
dotted and solid lines correspond to the symmetrical and the symmetry broken phase, respectively. The massless case is shown 
by dashed line. One can distinguish four different scaling regimes, namely A: U.V., B: symmetrical massive, C: precursor of 
condensation and D: spinodial phase separation regime. 

We identify at this point four different scaling regimes: 

A: In the U.V. scaling regime /c^ » 92(k): Gn ~ gn/k"^, the mass term is negligible. 

B: The explicit scale dependence disappears in the IR regime fc^ <C .92 (fc) of the symmetrical phase, where G„ « 
9n/ 92- This is a trivial scaling, there are no interactive relevant operators. 

C: The onset of the condensation is at G^^(fc^) = fc^ + g2{k) w and the higher loop effects start to dominate. 

D: The spinodal unstable region where the homogeneous ground state becomes unstable against infinitesimal fiuc- 
tuations, fc^ -|- 52 (fc) < 0. The leading scaling laws are coming from tree-level, classical physics. 

There are two crossover regimes. One is between the scaling regimes A and at A;^ ~ 52 (^)- Another one is 
between A and C when |$| < $infl- This will be extended to the region |$| < $vac by tree-level contributions. We 
shall look into the crossover of the symmetry broken phase. 

The best is to follow the evolution of the local potential the plane as shown in Fig. 8. The potential 

C/fe("S) should be imagined as a surface above the plane ($, fc^). The RG equations for the dimensionless coupling 
constants, g„ = /3„ where /3„ is given by Eq. (25) are integrated in the direction of the arrows. The initial condition 
Ua{4>) = m^(/)^/2fc^ + gk'^^'^cj)'*' set at A: = A, along the horizontal line on the top. As the running cut-off k is 
decreased and the RG equations are integrated the potential becomes known along horizontal line intercepting the 
ordinate at fc^. The spontaneous breaking of the symmetry in the vacuum implies that for sufficiently small 1$! 
the propagator explodes by approaching a singularity at A: = fccr(*I') as k decreases. This happens first at $ = 
when the running horizontal line touches the two curves shown in the Figure. These curves correspond to the 
generalization of the minimum and the infiection point of the potential for k^ > and are given by the relations 
k'^ + m% + g^l^iJQ + 0{h) =0 and k'^+m% + g^f^fi/2 + 0{h) = 0, respectively. We shall argue in section III C that 
the singularity lies close to the extension of the minimum. 
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fccr(^min(fc)) « k. 



(147) 



The numerical integration of the WH equation requires a truncation of the summation in Eq. (20) for the local 
potential. By using rimax < 22 the singular line can be located with reasonable accuracy. Close to the singular line 
on the plane ($, k'^) the coupling constants increase rapidly and the truncation of the potential is not acceptable. In 
the U.V. scaling regime (denoted by A in Fig. 7) the propagator is small, G = O (fc^), and the first term on the right 
hand side of Eq. (23) is dominant. As we approach k„ the propagator starts to increase, all term contribute equally 
in the beta functions and we enter into a new scaling regime (C in Fig. 7). 




<l> 



FIG. 8. The (<I>,/c^) plane where the local potential Uk{^) is defined. The potential should be imagined as a surface above 
the this plan. The initial condition gives C/a($) on the horizontal dashed line and the RG equation is integrated from this line 
in the direction of the arrows. The extension of the minimum and the inflection point of the potential from fc^ = to fc^ > 
is shown by sohd and dotted fines, respectively. 

In order to connect the two scaling regimes we shall consider fm,n{k/A) = df3m,k{4')/dgA,n, the dependence of the 
beta functions at k on the microscopical initial condition, imposed at fc = A. The simplest way to estimate fm,n{x) 
is to start with /m,n(l) ~ 1 which can be inferred from the form (23). The scaling law (145) allows us to write 

, (k\_ d~0M _ dPm{4>) (k\-^'-^ ,Jfk\-^'A 

up to corrections O (h), confirming the universality of the U.V. scaling regime. In fact, initial, microscopical value of 
the the irrelevant, non-renormalizable coupling constants {[g] < 0) leave vanishing trace on the dynamics as A/k oo. 

The numerical determination of the function fm,n{x) in d = 4 supports the prediction (148) in the U.V. scal- 
ing regime, namely the impact of a non-renormalizable coupling constants gn{A.) on the beta functions weakens as 
(fc/A)"[^"l for a; w 1. But this trend changes for .t < 1 as the critical value of the cut-off, fccr, found to be in agreement 
with Eq. (147), is approached in what / starts to increase. The potential must not be truncated close to the singular- 
ity and the numerical results appear inconclusive. But the increase of fm,n{x) at the crossover actually starts already 
far enough from k^ when all couplings are small, g„ << 1 and the truncation is safe. The lesson of this numerical 
result is that universality as used around the U.V. scaling regime alone is not valid anymore. The instability at the 
onset of the condensation at A; « k^ introduces divergences which are strong enough to overwrite the U.V. scaling 
laws and generate new relevant operators. Such an instability of the RG flow enhances the sensitivity of the physics 
at finite length scales on the microscopical parameters. 

The need of generalizing universality in a global manner has been shown for |$| < $vac- It remains to see if such a 
phenomenon can be observed at the highly singular point |$| = $vac, in the true vacuum. 

A similar phenomenon has already been noticed in connection with the BCS ground state. It has been pointed 
out that kinematical factors turn the four electron operator which is irrelevant according to the power counting into 
marginal for processes close to the Fermi level [60], [61] [62]. The coUinear divergences at the Fermi level drive the 
instability of the non-condensed vacuum and generate new scaling laws. Another similar mechanism is the origin 
of the strong long range correlations in the vacuum state of Yang-Mills models. They appear as a consequence of 
non-renormalizable, U.V. irrelevant Haar-measure term of the path integral [63] . 

We note finally that the view of the RG flow as a parallel transport, mentioned in section 11 C is particularly well 
suited for the studies of crossovers. This is because the sensitivity matrix (78) is a global quantity displaying clearly 
the sensitivity of the IR physics on the U.V. parameters by construction. 
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3. RG microscope 



We now embark on a rather speculative subject, playing with the possibility of matching different scaling laws. Let 
us start with the generic case, a model with U.V. and IR scaling regimes, separated by a crossover, as sketched in 
Fig. 6. We have a classification of the very same operator algebra at both scaling regimes. One can write any local 
operator ^ as a linear superposition of scaling operators, A = J2n''-A,nOn, the latter being the eigen-operators of 
the linearized blocking relations of a scaling regime, C'„(fc) w (A/fc)5'jC'(A)„. Let us denote by i/j^ the largest scaling 
dimension of scaling operators which occur in the a linear decomposition of A, va = max„z/„. We shall simplify 
matters by calling a local A operator relevant, marginal or irrelevant if i/^ > 0, = or < 0, respectively. We 
ignore marginal case by assuming that radiative corrections always generate non- vanishing scaling dimensions. 

One can distinguish between the following four cases: 

{^) ^u.vTiR^ relevant in both regimes, 

(b) ru.v.iiR- relevant in the U.V., irrelevant in the IR, 

(c) ixs.YXiK'- irrelevant in the U.V., relevant in the IR, 

(d) iv.v.im:- irrelevant in both regimes. 

The qualitative dependence of the corresponding coupling constants is shown in Figs. 9. 






(a) 



(b) 




(c) (d) 
FIG. 9. The qualitative dependence of the four types of coupling constants on the observational length scale. The U.V. and 
the IR scaling regimes are < a; < a;"" and a;"" < x, respectively. The dashed line in (c) corresponds to a different initial 
condition for the RG flow. 



Let us consider for example QED containing electrons, muons and ions with metallic density. There is an U.V. 
scaling regime for energies well above the electron mass. Suppose that the ground state is a BCS type superconductor 
where the IR scaling laws are driven by long range phonon mediated interactions. The electron mass is an example 
of the class (a). The muon mass belongs to class (b) owing to the decoupling of muons at low energy. By assuming 
that radiative corrections turn the four fermion vertex relevant it serves as an example of class (c). Finally the six 
fermion vertex is doubly irrelevant and belongs to class (d). 

The interesting class is (c). Suppose that we modify the values of the bare coupling constants at the initial condition 
and find the dashed line. Due to the universality of the U.V. scaling regime the two renormalized trajectories converge 
towards each other for < a; < x'^'" and run very close between the two scaling regimes, x^^ < x < x^^ . But in the 
IR scaling regime where the coupling constant is relevant the two trajectories start to diverge from each other. The 
question is whether the difference of the trajectories in the U.V. scaling regime has any impact on the difference found 
at the IR. It may happen that the IR phenomena is independent of the U.V. regime and the crossover smear completely 
out the small effects of the U.V. initial condition of the trajectory. But one can imagine that the extremely small 
differences left at the crossover remain important and lead to an initial condition dependent divergence at the IR. 
The answer is model dependent and can be given by the detailed analysis of the set of coupled differential equations. 
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In order to make this point clearer let us imagine a model with two parameters, say a mass m and a coupling 
constant g, both expressed in a dimensionless manner by means of the U.V. cut-ofF, the lattice spacing a. Suppose 
the following rather simple scaling laws: The beta function of the mass is approximately constant, 

CL dTTt 

7(m,3) = -— R^/3„>0 (149) 
m oa 



and of the coupling constant g can be written as 



dg 

P{m,g) = a— = 



m \ ( nicr \ 

X \ t^u.v. + X [ '^I.R. 

rricr J \ m J 



g (150) 



m(a) = m(ao) ( — ) , (151) 



where xi^) is interpolating smoothly between x(0) = 1 and x(oo) = 0. The running mass is given by 

and the length scale of the crossover between the U.V. and the I.R. scaling regimes is at acr ~ au.v.{''ncr / ''nu.v)^^ ^"^ ■ 
The asymptotical scaling for the coupling constant is 

/ \ —fu.v. 

g{au.v ) ~ cu.v.g{acr) ( ) > au.v.<. acr, 

g{ai.R.) « ci.R.g{acr) [^i^J , ai.R. > oicr- (152) 

Let us furthermore assume that I'u.v. < and i^lr. > 0, i.e. this coupHng constant belongs to the class (c). The 
sensitivity of the coupling constant in the asymptotical regions, g{au.v.) and g{ai,R,) on g{acr) is 

dg{au.v.) ( a,u.v. 

~ cu.v. 



dg{acr) 
therefore 

dgjai.R.) _ ci.R. ai'.i:a7.v: 

dg{au.v.) ~ cu.v. a^; « +"^-^- 



dg{acr) ' ' \ ttc, 

^^^"^■^■^~c,.«.f^y^" (153) 



(154) 



If the IR scaling regime is long enough then the difference between renormalized trajectories with different initial 
conditions imposed in the U.V. regime can be as large or even larger than at the initial condition. When the U.V. 
and the I.R. observational scales are related by 

/ ^ \ -VI.R./vu.V. 

au.v. ~ acr ( — ^ ) • (155) 

\ai.R.J 

then (154) is O (a°) ci.r./cu,v. and despite the 'focusing' of the universality in the U.V. regime the 'divergence' of 
the trajectories in the I.R. can amplify the extremely weak dependence on the initial condition at finite scales. The 
increase of ai,R, in this case enhances the sensitivity on the initial condition and we can 'see' the initial value of the 
non-renormalizable coupling constants at smaller distance a. 

The lesson of the numerical results mentioned above is that the scalar theory in the unstable regime realizes such a 
RG 'microscope'. The loss of the U.V. based universality as we approach the critical line on the plane ($, k'^) suggests 
that there is at least one relevant operator at the new scaling laws which is irrelevant in the U.V.. We do not know 
the eigen-operators of the linearized blocking relations in this region but this new relevant operator must contain the 
local monomials with n > 4. To make it more difficult, it is not obvious that this is a local operator as opposed 
to the scaling operators of an U.V. fixed point. 

In order to parametrize the physics of the scalar model we have to use this new coupling constant at the U.V. 
cut-off, as an additional free parameter. This is rather un-practical not only because we do not know the operator in 
question but mainly due to the smallness of this coupling constants in the U.V. and the crossover region. It seems 
more reasonable to use a mixed parametrization, consisting of the renormalizable coupling constants at the U.V. 
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cut-off and the new coupling constant taken close to the singularity where it has large enough value, ie fc w fccr- In 
the scalar model fccr ^ as the external source is turned off. We shall call the new coupling constants appearing in 
this manner hidden parameters. 

What operator corresponds to the hidden parameter of the scalar model? It seems reasonable to assume that the 
only observables whose value is determined at the onset of the condensate is just the magnitude of the condensate. 
Therefore the conjecture is that the strength of the c;ondensate, ^^oc; is dynamically independent of the renormalizablo 
parameters 52, 54 in d = 4. We have therefore the following possibilities in parametrizing the scalar model with 
spontaneously broken symmetry: 

• We identify a non-renormalizable coupling constant, gnr, which influences the value of the condensate and use 
the bare values of g2{^),gi{^),gnr{^) at the U.V cut-off. 

• One uses the renormalizable bare coupling constants from the U.V. end and the strength of the condensate from 
the I.R. regime, S'2(A), 54(A), ^^c- 

• One may use I.R. quantities only, say 52(0), 54(0), ^vac- 

Our point is that the scalar model with spontaneously broken symmetry has three free parameters instead of two. 
This rather surprising conjecture can be neither supported nor excluded by perturbation expansion because one can 
not connect the U.V. and the I.R. regimes in a reliable manner. The numerical studies of the scalar model carried out 
so far are inconclusive, as well, because they were constrained to the quartic potential and the possible importance of 
higher order vertices in forming the condensate was not considered. 

The modification of the scaling laws of the scalar model, necessary to generate a hidden parameter comes from the 
onset of a condensation. As discussed below in Section III C there are non-trivial tree-level scaling laws in a condensate 
which overwrite the loop-generated beta functions and may provide the dynamical origin of the hidden parameter. In 
a similar manner one can speculate about the role of the Cooper-pair condensate in the BCS vacumn. If there turns 
out to be a hidden parameter as in the one-component scalar model then the non-renormalizable effective coupling 
constants of the Standard Model would influence the supercurrent density in ordinary metals! Another possible 
example is the Higgs sector of the Standard Model, where the hidden parameter, if exists, would be an additional free 
parameter of the model. 

One may go further and inquire if dynamical symmetry breaking can modify the scaling law in a similar manner. 

The bound state formation, responsible to the generation of large anomalous dimensions in strong coupling QED [64] , 
in strong extended technicolor scenario [65,66], in the Nambu-Jona-Lasinio model [67] or in the top-quark condensate 
mechanism [68,66] may be a source of hidden parameters, as well. 

Hidden parameters represent an unexpected 'coupling' between phenomena with very different scales and question 
of our traditional strategy to understand a complex system by analyzing its constituents first. The axial [70] or scale 
[71] anomalies represent a well know problem of this sort, except that the regulator independence of the anomaly 
suggests that for each classical IR fixed point there is a single 'anomalous' one without continuous fine tuning. 

Universality is expressed in the context of High Energy Physics by the decoupling theorem [69]. Let us start with 
a renormalizable model containing a heavy and a light particle and consider the effective theory for the light particle 
obtained by eliminating the heavy one. There are two possible classification schemes for the effective vertices for 
the light particle, generated by the elimination process. One is according to the U.V. scaling laws in the effective 
theory, i.e. there are renormalizable and non-renormalizable vertices. Another scheme is based on the strength of 
the effective interactions, i.e. there are vertices which are stay constant or diverge when the ratio of the light and 
the heavy particle mass tends to zero and there are vertices whose coefficients tend to zero in the same limit. The 
decoupling theorem asserts that these two classification schemes are equivalent, i.e. the non-renormalizable effective 
coupling constants are vanishing when the heavy mass diverges. The existence of the hidden parameter is, at the final 
count is a violation of this theorem. 



4- The Theory of Everything 

The appearance of a hidden coupling constant renders our goal of understanding the hierarchy of interactions in 
Nature extremely complicated. In fact, the RG flow of the Theory of Everything must be imagined in the space of 
all coupling constants in Physics, including the parameters of the Grand Unified Models down to different models 
in Solid State Physics or hydrodynamics in the classical regime, as indicated in Fig. 10. Different fixed points are 
approached by the trajectory when the scale dependence is dominated by a single interaction. But all fixed point 
except the U.V. one is avoided due to the non-renormalizable, irrelevant vertices [69] generated by the dynamics of 
the next fixed point in the U.V. direction. 
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Such a wandering in the space of coupling constants echoes the age old disagreement between High Energy and 
Solid State Physics. It is usually taken for granted in the High Energy Physics community that the sufficiently precise 
determination of the microscopical parameters of the Theory of Everjrthing would 'fix' the physics at lower energies. 
The obvious difficulties to extract the high energy parameters from experiments which render each new experiment 
expensive, long and large scale operation indicate that this might not be a practical direction to follow. In other 
words, the relevant coupling constants have positive Lyapunov exponent and render the trajectory extremely sensitive 
on the initial conditions. Therefore the characterization of the RG flow of The Theory of Everything by the initial 
condition, though being possible mathematically, is not practical due to the finite resolution of the measurements. 
The other side of the coin, the physical parameters of a fixed point are clearly the relevant coupling constants of the 
given scaling regime. They are the object of Solid State Physics, as far as the scaling regimes QED, CM and IR are 
concerned in Fig. 10. One looses sight of the fundamental laws but gains predictive power by restricting oneself to 
local studies of the RG flow. 




FIG. 10. The rcnormalization group flow of the Theory of Everything. The branching drawn at an IR scale is a phase 
transition driven by the environment, such as particle or heat reservoirs. 

C. Instability induced renormalization 

In the traditional applications of the RG method which are based on the perturbation expansion the RG trajectory 
is sought in the vicinity of the Gaussian fixed point. One can describe in this manner crossovers which connect scaling 

regimes which share the same small parameter. We can enlarge the classes of accessible crossovers by relying on the 
saddle point expansion, the only systematical non-perturbative approximation scheme. The generic way to induce 
new crossover in this scheme is by passing a condensation at a finite length scale [72] . 

One may object that there is nothing surprising or new in finding that the physics is fundamentally changed by 
a condensation mechanism. What additional knowledge can then be gained by looking into condensation by the 
method of the RG? The point is that the classification of operators around a fixed point which was achieved in 
the framework of the perturbation expansion changes essentially by the condensation, ie the appearance of a non- 
trivial saddle point. To sec how this happens we return to the strategy followed at Eq. (2). Let us consider an 
observable A in a model with a single coupling constant g for simplicity, computed in the saddle point expansion, 
{A) = Fo (5, A) + hFi (5, A) + O (Jt?) , where [g, A) denotes the £-loop contribution. We can obtain the beta function 
by taking the derivative of this equation with respect to cut-off A, 

The loop corrections F£{g,A), £ > 1 are polynomials of the coupling constants g. But the leading order, tree-level 



P = kdkg = - 



kdk{Fo + nFi) 
dJFo + hFi) 



+ o (n^) = - 
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piece usually has stronger dependence on g and may become singular as g 0. The tree-level contribution may induce 
qualitatively new scaling laws with new set of relevant operators. Due to the singularity at g = the condensation 
actually realizes the dangerous irrelevant variable scenario [73]. 

At the onset of an instability, region C in Fig. 7 certain modes experience strong nonlinearity and develop large 
amplitude fluctuations. But as we enter in the unstable regime one may hope to recover quasiparticles with weak 
residual interaction after having settled the dynamics of the unstable modes. At least this is what happens within the 
framework of the semiclassical, or loop expansion where the fluctuations remain small after the proper condensate is 
found. There are models where the vacuum state is just at the edge of the instability, like the scalar model discussed 
above in the symmetry broken phase. 

The tree-level contributions to the blocking transformation, if exist, are more important then the loop corrections. 
This is the reason that we shall consider the tree-level, O (fi ''^ RG flow. There is no evolution in this order so long the 
system is stable and the saddle point is trivial. But as soon as we arrive at the unstable region the RG flow becomes 
non-trivial. 



1. Unstable effective potential 

Before turning to the actual tree-level blocking let us review briefly the kind of instabilities one expects. As we 
enter into the unstable regime k'^ + 52(A) < the trivial saddle point, 0^ = becomes unstable in the blocking (7). 
To understand this instability better we introduce the constrained functional integral, 

Z^ = J l?[</.]e-^-[^l5 {vj/^' ^^^^^ 

corresponding to a conserved order parameter. By following the RG trajectory until the IR end point one eliminates 
all but the homogeneous mode by keeping the partition function unchanged, 

Z$ = e-*^=°f*i = e-^^'==°(*^ (158) 

To find another local representation for this constrained partition function we consider the generator functional 
density w\j] = W[j]/V defined by means of Eq. (5) and compute its Legendre transform, Vqs{^) in the mean field 
approximation. For this end we set jx = J and write 

-14ff($) = w[J] - J$ = mm[-Uk=oi<l>) + J{<t> - *)]■ (159) 

The minimization with respect (j) and J yields J = dw{J)/dJ and and = <&, respectively. Finally we have 
Ves{^) = Uk=o{^), the result announced after Eq. (16) so long the mean field approximation is reliable, namely 
in the thermodynamical limit and in the absence of large amplitude fluctuations. It is further known that Veff{^) is 
convex. 

The instabilities of the kind mentioned above are well known in the case of first order phase transitions. But similar 
instabilities may appear at higher order phase transitions, as well. In fact, the magnetization of the Ising model as 
the function of the external magnetic field shows discontinuous behavior below the critical temperature. Therefore 
the free energy constrained into a sector with a given magnetization displays such kind of instabilities. 

Suppose that the Uk=o{^) determined perturbatively has degenerate minima and a concave part. Then it is 
advantageous to introduce two curves on the ($, fc) plane of Fig. 8, $infl(fc) and $min(fc) defined by fc^$min + 
f^fe^^(^min) = 0, and fc^ -t- U^\^in{i) = 0. The stable region in the mean-field approximation is |$| > i>inin(0) = 
$vac- For $infl(0) < |$| < $min(0) there are two minima in (p for the last equation in (159). One of them is 
metastablc, ie is unstable against sufficiently large amplitude modes. The spinodal phase separation, the instability 
against infinitesimally small amplitude fluctuations occurs when some eigenvalues of the second functional derivative 
of the action becomes negative. This is region D in Fig. 7, |$| < $infi(0). It is unreachable by the mean- field 
treatment because no local minimum is found in (j) which would satisfy \(f>\ < $infi(0). To understand better the 
nature of these unstable regions one has to go beyond the mean-field approximation and to follow the dynamics of 
the growing inhomogeneous instabilities. We shall use the tree-level WH equation to deal with such large amplitude, 
inhomogeneous fiuctuations. 
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2. Tree-level WH equation 



The construction of the condensate within the unstable region impHes the minimization of the action with respect 
to a large number of modes. The RG strategy offers an approximation for such a rather involved problem, it deals 
with the modes one- by-one during the minimization [75], [74]. 

The negative curvature of the potential makes the saddle point in (9) non-trivial. 



[<^] = mjn Sk [</. + 0'] =Sk[(t> + ki [<!>]] + Sk [</.] . (160) 



Since the saddle point depends in general on the background field, (f)ci = <pci[(t>]j the action is modified during the 
blocking and we find a non- trivial RG flow [74]. The local potential approximation to the tree- level blocking is 



Uk-Ak{^) = min / 

4>cl.x J X 



(161) 



One encounters here a conceptual problem. The fc-dependence of the trcc-lcvcl RG flow might well be singular 
since there is no obvious reason that the saddle points which are usually rather singular functions of the parameters 
evolve smoothly. This problem has already been noticed as a possible 'first order phase transition' in the blocking 
which induces discontinuous RG flow [76], [77], [78]. It was later shown by rigorous methods that the RG flow is 
either continuous or the blocked action is non-local [79]. The resolution of this apparent paradox is that the saddle 
point actually develops in a continuous manner as we shall shown below. As the RG flow approaches the onset of 
the condensate, k kcr then the new scaling laws generate such an action that the saddle point turns out to be 
continuous in k. 



3. Plane wave saddle points 

The saddle points, the minima in (160) satisfy the highly non-linear Euler-Lagrange equations whose solutions are 
difficult to find. The use of sharp cut-off slightly simplifies the problem since it reduces the functional space in which 
the minimum is sought to <f) & J^k\^k-Ak- We shall retain the the plane wave saddle points only, 

(t>ci,p = ^ [e*^'°'5p,fcefc + e^^''''6p-kek] ) 4'ci,x = Pk cos(fcefc • x + Ok)- (162) 

The parameter Ok and the unit vector efe correspond to zero modes, they control the breakdown of translational and 
rotational symmetries. 

This is a key to what happens at the unstable line: Despite the discreteness of the internal symmetry (j) —(j) there 
are soft modes because the inhomogeneous saddle points break the continuous external symmetries. The condensation 
of particles with non- vanishing momenta automatically generates Goldstone modes. This phenomenon is well known 
for solids where the saddle point is a crystal of solitons which breaks external symmetries and there is no integration 
over the zero mode to restore the symmetrical ground state. Our point is that the soft modes make their appearance 
even if the symmetry of the ground state is restored by the integration over the zero modes. 

The amplitude pk is determined by minimization, 

C/fc-Afe($) = min ( ifc2p2 + i /" dyUk{^ + Pkcosy)] , (163) 

in the local potential approximation. The numerical implementation of this iterative procedure to find the RG flow 
with an initial condition imposed on the potential at fc = A gave the following results [74]: 

1. The recursive blocking relation, Eq. (163) is not a finite difference equation. Despite of this the flow converges 
as Afc 0. 

2. The amplitude of the saddle point satisfies the equation $ -|- p{k) = i>niin(A:). 

3. The potential obtained by the tree-level blocking is Uk{^) — —^k"^^"^. 

4. The tree-level results above hold independently of the choice of the potential at the cut-off. 
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The key is point 3, point 1 follows immediately from it. The lesson of this result is that the action is degenerate for 
the modes at the cut-off, the kinetic and potential energies cancel. The 'best' effective theory for a given plane wave 
mode is the one whose cut-off k is slightly above the wave vector of the mode. This result suggests the degeneracy of 
the action within the whole unstable region. Conversely, if we can show that the action is degenerate at the cut-off 
within the unstable regime we established this potential. 

We prove by induction in the number of steps k k—Ak that the variation of the action density within the unstable 
regime is O {Ak). This result protects the consistency of the saddle point expansion for d > 1 since Afc >2'k/L where 
L is the size of the system and therefore the variation of the action is O (i'^^^). Recall that kcr denotes the cut-off 
where the kinetic snd the potential energies cancel. 

First step: Let us denote by k' the value of the cut-off at the first occurrence of non-trivial saddle point in the 
numerical implementation of Eq. (163). It obviously satisfies the inequality k^ — Ak < k' < k^r- The trivial 
saddle point (pes — becomes unstable for blockings with |$| < $infl(fc') in which case |$ + 0cs,a;| < ^inf[{k')^ ie 
\4>cs,x\ = O (^^/Akj . The term O (j)^ is canceled in the action on a homogeneous background field so the O (^cs) 
contribution gives O {Ak) variation. 

Induction: Suppose that the variation of the action is 0{Ak) in the unstable region and we lower the cut-off, 
k ^ k ~ Ak. At the new cut-off the balance between the kinetic and the potential energy is lost by an amount of 
O {Ak) since the potential energy is still the given by k but the kinetic energy corresponds to the lowered cut-off, 
k — Ak. Thus the negative curvature potential energy wins and the action bends downward in the unstable region as 
the function of the amplitude of the plane wave. By assuming that the amplitude is stabilized at an O (Afc*^) value 
(will be checked later) the variation of the action density O {Ak) in the mistablc region. 

Point 2 can be understood in the following manner: During the minimization of the action pk slides down on the 
O {Ak) slope until the potential starts to increase again. We can find where this happens by equating the slope 
of the original, bare potential with that stated in point 3. One would expect that the deepest point is reached at 
Pk = *I'infi(fc)- But the O {(f)'^) term in Uk{(j)) is not yet strong enough to make the potential increase strong enough 
here. The two slopes agree just at pk = $min(A;). 

Result 4 which follows from 3, as well, can be called super universality since it reflects scaling laws where all coupling 
constants are irrelevant. Note that the fc-dependence is continuous through the whole RG trajectory. 

Point 3 is a generalized Maxwell construction. It reduces to the traditional Maxwell construction for fc = 0, to 
the degeneracy of Uk=o{^) for |$| < $rnin(0) = ^vac- The naive Maxwell cut, applied for the concave part of the 
potential would give the degeneracy for |<I>| < <I>infi only. The problem with this argument is that it produces an 
effective potential which is convex everywhere except at |(f>| = $infi(0). The second derivative of the potential is ill 
defined and the first derivative is discontinuous at this point. By placing the system into a finite box the singularity 
is rounded off and the second derivative becomes finite, but turns out to be negative in a vicinity of |$| = $infl(0). 
Convexity regained only if the cut is extended between the minima, |i>| < $min(0). 

The generalized Maxwell construction in the mixed phase of a first order phase transition can be understood by 
the dynamics of the domain walls. The fiatness of certain thermodynamical potentials in the mixed phase reflects 
the presence of zero-modes, the location of domain walls. Such a rather simple kinematical mechanism which is 
independent of microscopic details is the source of the 'super universality', point 4 above. The role of the domain 
walls is played by saddle points in our computation, the cosine function in (162) realizes infinitely many equally 
spaced, parallel domain walls. The integration over the zero modes 6k and according to the rules of the saddle 
point expansion restores translational and rotational symmetries of the ground state and reprodiicc^s the mixc^d phase. 

One might object our independent treatment of the plane wave saddle points. The RG equation (163) handles the 
plane waves in the consecutive momentum space shells independently which seems as a vast oversimplification. But 
one should recall at this point that according to the general strategy of the RG method the dynamics of the modes 
eliminated during the blocking is retained by the modification of the effective coupling constants. This is not always 
possible since there arc more modes than coupling constants, the problem mentioned at the end of section II A 1. 
More precisely, the general framework of the blocking what is in principle always applicable is turned into a powerful 
scheme when an approximation is made. One truncates the effective action and assumes that the solution of this over- 
determined problem exists. This is the main, and so far unproven, assumption of the RG method. Accepting this point 
we can determine the blocked action by means of a homogeneous background field and use it for non-homogeneous 
field configurations at the next blocking step as written in Eq. (163). 

4. Correlation functions 

One can compute the tree-level contributions to the correlation functions in a combination of the mean-field and 
the saddle point approximation [80]. Let us split the complete functional integral into the sum of contributions which 
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come from the unstable and the stable regions. The fluctuations from the stable regions are governed by a non-trivial 
action and they are taken into account as loop corrections. The fluctuations in the unstable region experience a flat 
action and their contributions must be taken into account on the tree-level. The complication of integrating up these 
contributions is the determination of the region where the action is flat. Our approximation consists of estimating 
this region for each plane wave independently, ie extending the integration over the amplitude Vp of the degenerate 
plane wave 



4>ci{x) = rp cos{pepX + 0p) 



(164) 



over the interval —pp < r < pp, where pp is given by (162). Let us denote the integration over the resulting domain 
by [r] and write the correlation function in momentum space in this single- mode determined flat region as 



J D[0]D[e]D4r] J D[e]D[e]D4r]rprg 



(165) 



Due to the flatness of the action these integrals are purely kinematical and can easily be carried out. All integration 
whose variable does not show up in the integrand drops out. The integration over the shift of the plane waves, the 
phase angle 0, restores the translation invariance. 
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where $min(^) defined as $min($min(^)) = ^- The Fourier transform of this propagator describes a diffraction 
type oscillation with characteristic length scale ^ = k~^, the characteristic feature of domains in a homogeneous state. 



5. Condensation as crossover 



The scaling laws of the scalar model change already approaching the condensation in the stable region C of Fig. 7. 
Inside the instable region D the action with the potential given by point 3 of section III C 3 becomes invariant under 
blocking, ie the whole region D is an IR fixed point. 

The vacuum of the scalar model is a single, homogeneous coherent state consisting of zero momentum particles 
when |$| > $vac- When l^l < ^vac we encounter the singularities as seen above. This singularity may not be real, it 
might be smoothened out by higher order vertices and a saddle points appearing for fc < fccr- But the vacuum is the 
superposition of inhomogeneous coherent states with characteristic scales £ > 1/fccr- Each of them breaks space-time 
symmetries but their sum remains symmetrical. 

One may wonder if any of the discussion applies to the true vacuum |$| = ^vac- For a weakly coupled system 
the fluctuations are small and one may hope that the fluctuations around $vac arc stable and un-influenced by what 
is happening with the unstable modes < ^vac- But the answer to this question is rather uncertain and is open 
because the true vacuum is just at the border of the instabilities and the typical fluctuations around the true vacuum 
{(j>x — ^vac) ~ 5'^(0)/5'4(0) penetrate into the unstable region. 

The issue of the modification of the scale dependence by a condensate [64] is rather general and points far beyond 
the simple scalar model. 

The vacuum of asymptotically free models supports strong correlations at long distance. In the case of 4 dimensional 
Yang-Mills models the naive, perturbative vacuum is destabilized by the one-loop level effective action which predicts 
that the vacuum is made up by a coherent state of gluons [82]. But this state cannot be the true vacuum. The problem 
is not only that the field generated by the condensate is strong and spoils the saddle point expansion but it turns out 
that there is an unstable mode. The true vacuum is supposed to be found at even lower energy densities where the long 
range fluctuations restore the the external and color symmetry broken by a homogeneous condensate of the charged 
spin one gluons [83] and the vacuum is thought to contain domains of homogeneous field in a stochastic manner [84]. 
This scenario is close to the view of the mixed phase of the scalar model developed here with the difference that the 
instability comes from the loop or the tree-level renormalized action in the Yang-Mills or scalar model, respectively. 
Furthermore the instability in the Yang-Mills model can be avoided by extreme environment only, by immersing the 
system into strong external field or bringing into contact with heat or particle reservoir. 
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Similar instability occurs in vacua containing condensate of bound states. The BCS vacuum is made homogeneous 
in a non-trivial manner when represented in terms of the electrons making up the Cooper pairs. The spontaneous 
breakdown of the chiral invariance in QCD manifests itself in the condensate of quark-anti quark pairs which is 
homogeneous after integrating out the instanton zero modes only [85]. The crossover from the U.V. scaling regime to 
the instability takes place at few l/^ where £ is the size of the bound states. Higher order derivative terms appearing 
as effective vertices may generate similar crossover, as well [72], [86]. 

Finally it is worthwhile mentioning the tunneling phenomena, the dynamical extension of the instabilities considered 
above. The interesting feature of the dynamical realization of the tunneling by means of time-dependent, tree- level 
instabilities is that they take place in conservative systems but the long time distribution agrees with the equilibrium 
predictions coming from the canonical ensemble [81]. 

D. Sine-Gordon model 

It is not unusual that one arrives at periodic or variables in the construction of effective theories, cf non-linear 
sigma models or non-Abelian gauge theories. Such variables represent a double challenge. One problem which can be 
considered as local in the internal space is that the perturbation expansion around a minima of the potential should 
keep infinitely many vertices in order to preserve the periodicity. Another, more difficult, global problem is a conflict 
between the two requirements for the effective potential for periodic variables, namely periodicity and convexity. One 
can fulfill both requirements in a trivial manner only, by constant effective potential. If true, this conclusion has far 
reaching consequence for the phenomenological description of such systems. As a case study we shall consider the 
simplest model with periodic variable, the two dimensional sine-Gordon model [87] , [88] . 

1. Zoology of the sine- Gordon model 

The sine-Gordon model 

LsG = \{d^cl)f + ucos(ici> (167) 

has been shown to be equivalent with the X-Y model [89], [90], [87], the Thirring model [91], [92] and a Coulomb 
gas [93]. The methods of Refs. [90], [91] and [92], approximate duality transformation, bosonization and semiclassical 
approximation are valid in certain regions of the coupling constant space. The maps used in Refs. [87] and [93] are 
exact. The mapping [87] Vx = e.^'^" transforms the model (167) into the compactified sine-Gordon model 

icsG = ^a^V*5^^ + |(V' + V'*)- (168) 

The models (167) and (168) are equivalent in any order of the perturbation expansion in continuous space-time where 

the configurations are assumed to be regular. 

The X-Y model which is characterized by the action 

-Inz^m^, (169) 

X 

where T = (fi , h/T = u, z \s the vortex fugacity and nix denotes the vortex density. The RG equation of the X-Y 
model obtained in the dilute vortex gas limit [89] , [90] 

a— ^ ir^z^ -nT^h^, a— = I 2 ] h, a— = 2 ] z. (170) 

da ' da \ 4nJ ' dz \ tJ ^ ' 

The X-Y model appears as a generalization of the compactified sine-Gordon model since the lattice regularization 
transforms the Lagrangian (168) into (169) with z = 1. What is surprising here is that the vortex fugacity is a relevant 
operator in the high temperature phase of the X-Y model but it is entirely missed by the lattice regularization! Does 
that mean that the sine-Gordon model Lagrangians (167) and (168) are incomplete? 

The vortex term could be missed by the regulator because the plane z = is RG invariant. Therefore it is consistent 
to exclude vortices from the path integral in the sine-Gordon model given in continuous space-time. When we recast 
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the model in lattice regularization without ever thinking about vortices our lattice action neither suppresses, nor 
enhances the singular configurations and sets z = 1. But once the vortices are not suppressed, 2: 7^ the renormalized 
trajectories moves away from from the plane of fixed z. There is still no problem in the low temperature phase because 
all what happens is that an irrelevant coupling constant became fixed. This is how regulators work. But the problem 
is more serious in the high temperature phase. The Lagrangians (167) or (168) can not give account of this phase 
because they lack a relevant, renormalizable coupling constant. The perturbative continuum theory without vortices 
is consistent but once vortices are available kinematically their density evolves in a non-trivial manner. 

There is no problem with the lattice regulated X-Y model because the vortex fugacity gets renormalized due to 
the cut-off dependence of the vortex action. In other words, we can move along the RG flow by adjusting T and h 
only, the renormalization of z is carried out 'automatically' in the partition function. The renormalization of 2; is a 
problem in the continuous formalism only where the vortices are introduced formally as point-like charges. 

This is an unexpected mechanism which brings the singular nature of the field configurations into play during 
renormalization and may plague any quantum field theory. The knowledge of the classical action in continuous space- 
time leaves open the possibility that we have of adjust the fugacity of certain localized singularities or topological 
defects not considered in the continuum. 

It is interesting to speculate about similar phenomenon in QCD. Lattice QCD is constructed with Inz = where 
z is the fugacity of some topological defects, like instantons, monopoles, merons etc. These topological defects and 
singularities are supposed to play an important role in the confining vacuum, as vortices do in the high temperature 
phase of the X-Y model. The renormalization of the fugacities of these objects should first be studied in lattice QCD 
in order to construct a continuum description of the vacuum. 

The sine-Gordon model possesses a topological current, j^^^ = /3e/ji/C^i/'/>a;/27r which is obviously conserved when the 
path integral is saturated by field configurations with analytic space-time dependence. Its flux, the vorticity, gives the 
soliton number. The world lines of the sine-Gordon solitons end at the X-Y model vortices, cf Fig. 11. The distance 
between the bound vortex-anti vortex pairs shrinks with the lattice spacing in the low temperature continuum limit. 
Any measurement with finite resolution loose sight of the instability of solitons in the renormalized theory. The 
average distance between vortices stays flnite, cut-off independent in the continuum limit of the high temperature 
phase and the soliton decay can be observed. One expects the breakdown or at least important modiflcation of the 
bosonization transformation in this phase. In fact, the non-conservation of the topological current requires fermion 
number non-conserving terms in the fermionic representation, a fundamental violation of the rules inferred from the 
weak coupling expansion. 



V 




FIG. 11. The world line of a soliton starts at a vortex (V) and ends at an anti- vortex (A) in the two dimensional space-time. 



2. Effective potential 

In order to understand better the dynamics of the long wavelength modes we compute the effective potential in 
the sine-Gordon model. We follow the WH method truncated to the local potential approximation which picks up 
the tree level evolution, too. It is obvious that period length of the potential 2n/(3 remains RG invariant in this 
approximation. The tree ands the loop level RG equations are 



Uk-Ak{<f>) = min 
p 



— -I- J duUk{(p + pcos{nu)) 



(171) 



and 
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in the plane z = 0. 

We use the Fourier expanded form for the potential 
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Uk{(t>) = Un{k) COS {nP(f)) 



(173) 



n=0 



and compute the leading order contribution to the WH equation when expanded in the potential, 



kdkU, 




(174) 



in terms of the the dimensionless coupling constants u„ = Un/k"^- This agrees with the second equation in (170). The 
solution of (174) is 



The coupling constants u„ are irrelevant in the disordered phase T = > Stt and the effective potential obtained 
for A; = is flat. The coupling constants n < Sw/T are relevant in the ordered phase and the effective potential is 
non-trivial. At this pointe one suspects that Maxwell construction interferes with the evolution because a non-trivial 
periodic functions necessarily has concave regions. To settle this question one has to rely on the numerical solution of 
the evolution equation (171). We followed the loops-generated evolution (172) from the initial condition set at /c = A 
until the stability is lost at A; = ka for background field values which lie at local maxima of the periodic potential. 
For k < fccr equation (171) was used. The result is that the coupling constants approach zero as A; ^ due to the 
plane wave saddle points. The effective potential is trivial in both phases. 



The local potential Uk{(j)) flattens without developing singularities in the high temperature disordered phase. The 

transformation 



is a discrete symmetry of the action and is preserved in the vacuum since the potential barrier between the minima 

is vanishing for long range modes. Since we are at the lower critical dimension, d = 2, the large amplitude long range 
modes realize the 'tunneling' between the minima of the periodic potential. On the contrary to this situation, the 
potential develops discontinuous second derivatives in the ordered, low temperature phase and the instability driven 
flattening of the potential Uk{4>) reflects the survival of barriers between the minima of the potential. All this looks 
like a spontaneous symmetry breaking, therefore our conclusion is that the transformation (176) is not a symmetry 
of the vacuum in the ordered, low temperature phase. 

The usual circumstance under which such a phenomenon arises is the multiple connectedness of the internal space. 
In the present context the non-linear U{1) cr-model parametrization, Eq. (168), is based on the internal space U{1), 
with the fundamental group Z generated by the transformation (176). The dynamical breakdown of the fundamental 
group symmetry is a genuine quantum effect. In fact, the path integral formally extends over all homotopy classes, 
this is the symmetrical phase. When the dynamics develops sufficiently high barriers between the homotopy classes 
the path integral becomes restricted to a single homotopy class. 

It is important to recall that the time evolution described by the Schrodinger equation can be derived from the path 
integral by performing infinitesimal variations on the end point of the trajectory. Therefore the consistency of the 
dynamics can be maintained when the path integration is restricted into any functional subspace which is closed under 
continuous deformation of the trajectories. The loss of the interference between homotopy classes is the characteristic 
feature of the symmetry broken phase. 

Such symmetry breaking is the key to understanding the way quarks become deconfincid at high temperature or 
the droplet phase is formed for quantum liquids [94]. The configuration space for global gauge transformations is 
SU{3)/Zt, for gluons, with the fundamental group Z3. The quark propagator vanishes in the symmetrical phase due 
to the destructive interference between the three homotopy classes. In the Z3 symmetry broken deconfined phase 
there is no interference and quarks can be observed. Quantum liquids in the first quantized formalism display similar 
symmetry breaking pattern. The coordinate space for iV particles is R^^ /Sn where Sn consists of permutations of the 
particles. The absence of the overlap among states belonging to different droplets suppresses the (anti)symmetrization 
of the states and reduces the exchange symmetry of the ground state. 




(175) 



3. Breakdown of the fundamental group symmetry 




(176) 
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How can wc recognize the dynamical breakdown of the fundamental group symmetry? The simplest strategy, to 
look for the minima of the effective potential fails because Maxwell construction hides any structure in a periodic 
potential. The answer to this question lies in the topological structure of excitations and the order parameter will be 
a topological susceptibility. 

This is rather natural since the phase transition is the restriction of the path integral into a single homotopy class, 
the imposition of a topological constraint. Lot us assume that the analysis sketched above for the plane z = remains 
qualitatively valid for z ^ 0,ie the effective potential is always flat, due to the loop-generated evolution or the Maxwell 
construction in the disordered or ordered phase, respectively. The topological invariant characterizing the homotopy 
classes is the soliton number. In the weak coupling expansion (ordered phase) the path integral is constrained into 
a single homotopy class hence the soliton number is conserved, its susceptibility is vanishing in the continuum limit. 
The stability of the soliton, based on the continuity of the time evolution is lost in the disorder phase because the short 
distance, large amplitude fluctuations of the typical field configuration extends the path integral all soliton number 
sector. Therefore the susceptibility of the soliton number is non- vanishing in this phase. 

4- Lower critical dimension 

It is worthwhile noting the manner the Mermin-Wagner-Coleman theorem [95] appears in the framework of the local 
potential approximation. By setting (72 = for d = 2 the linear part of P„ in Eqs. (23) provides a finite evolution but 
the higher order terms make the renormalized trajectory to diverge as fc ^ unless the coupling constants approach 
zero. We can construct interactive, massless models in two dimensions so long the running coupling constants approach 
zero sufficiently fast in the IR regime as in the case of the sine-Gordon model, presented above. The retaining of the 
terms O (9^) in the gradient expansion provides another mechanism to suppress the IR divergences by generating 
singular wave function renormalization constant, a,s k —^ [96]. 

In both cases the way to avoid the non-interactive system predicted by the Mermin-Wagner-Coleman theorem is to 
use the RG improved perturbation expansion. In the ordinary perturbation series we have the contributions like 

5" / fp^,-,Pk (177) 

Jpir-'iPk 

and the integral is IR divergent in rf = 2. One has running coupling constants, g gp, in the RG improved 
perturbation expansion where the analogous contribution is 

/ 9p(p,,...,p^)fpi,-,Pk- (178) 

•'Pl,----Pk 

When the IR fixed point is Gaussian, limp^o5P = 0, then the theory remains interactive and the running coupling 
constants suppress the IR divergences. 

E. Gauge models 

The blocking step of the RG scheme, being scale dependent, violates local symmetries which mix modes with 
different length scales. 

The most natural way to deal with this breakdown of gauge invariance is to cancel the gauge non-invariant terms 
by fine tuning non-invariant counterterms [97]. But this is possible in perturbation expansion only. Another, better 
suited strategy to the RG method is to argue that the way Ward identities are violated by blocking shows that the 
BRST invariance is broken by the running cut-off only and the BRST symmetry is recovered at the IR fixed point 
when the cut-off is removed [20] , [98] . This were certainly a good procedure if it could be implemented without further 
truncation. In a weakly coupled model such as QED one could construct approximations with controllable errors. 
In asymptotically free models like QCD non-perturbative long range correlations make any scheme which imposes 
gauge invariance in an approximate manner unreliable [63]. The demonstrate this point consider the static force law 
between two test charges in a Yang-Mills theory in a scheme where the gauge invariance. Gauss' law is implemented 
approximately. The gauge non-invariant components of the vacuum state appear charged, by definition. Thus we 
have an uncontrollable color charge distribution around the test charges. When the distance between the test charges 
are large enough then it will be energetically favorable for the uncontrollable charges to break the fiux tube. Wc loose 
the string tension in a manner similar how it happens in the true QCD vacuum due to the virtual quark-anti quark 
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polarizations. In more formal words, an arbitrarily weak gauge-dependent perturbation can change the long range 
features of the vacuum as in the ferromagnetic phase of the Ising model. 

The most natural way to guarantee gauge invariance is to achieve independence on the choice of gauge by using 
gauge invariant quantities only. The proposal of Ref. [99] goes along this line but the use of loop variables renders 
the computation rather involved. We discuss now a version of the RG scheme in the internal space which at least 
for Abelian models produces the effective action without any gauge fixing. The electron mass combined with chiral 
transformation has already been used in generating evolution equation in Ref. [100]. We shall follow a simpler and 
more general by avoiding gauge dependence in an explicit manner for the photon Green functions [101]. The gauge 
invariant electron composite operators are not difficult to include [41]. 



1. Evolution equation 



Let us consider the generator functional 



(179) 



where — df^ + lA^ and the dimensional regularization is used to render W finite. We control the quantum 
fluctuations by modifying the action, S ^ S + S\ with 



(180) 



where T^^, = g^^ — d^d^/U. The infinitesimal change of the control parameter A — > A + AA modifies the free photon 
propagator 



D^'^x -y)^ D'^'ix - y) + J dzDf'Pix - z)n,Tp^D''''{z - y) + O {{AXf) 

in a manner reminiscent of the Callan-Symanzik scheme. 

The evolution equation (108) with the present suppression takes the form 



(181) 



d,T^[A] = — tr 



5AJA^ 



(182) 



in terms of the photon effective action Tx[A] = T\[A] + Sx[A]. We project the evolution equation on the functional 
space given by the ansatz 



where 



f[A] = i ^ A^^,D-\id)T^''A,^, + C[A] 



(183) 



(184) 



_ _ rpixv Q j^j jg ^ gauge invariant functional. 
The control parameter A 'turns on' the fluctuations of the photon field. Therefore the electrons loop contributions 
to the effective action must already be present at the initial condition which is chosen to be C\A\ = —itTln{iI^ — m) 

at A = Aq. 

The second functional derivative matrix is written as 



mA] 



5Afj,5A^ 5Afj,5A^ 
and the inversion is carried out by expanding in the non-diagonal pieces to write the evolution equation as 



n=0 



(185) 



(186) 
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2. Gauge invariance 



We show now that the Umit a — > can be taken in the evolution equation without hitting any singularity. The 
gauge fixing parameter a enters through the photon propagator D in the Ncuman-cxpansion of the right hand side of 
Eq. (186). The a-dependent longitudinal contributions of the first and the last D factor are suppressed by the gauge 
invariance of the suppression term, represented by the transverse projection T'^" in (186). The longitudinal photon 
contributions of the internal propagators are suppressed by the gauge invariance of the effective action, 

a,^=0. (187) 

According to this equation 5^C [A\ / SA^^SA^ does not mix the longitudinal and the transverse modes, ie non- longitudinal 
contributions appear in the evolution equation. 

We choose a gauge in the argument above what was relaxed at the end of the computation. But the steps followed 
remain well defined even if we start with a = 0, without any gauge fixing. Our argument about the decoupling of the 
longitudinal and transverse contributions to the evolution equation still applies but it is not clear if the longitudinal 
part was well defined. This subtle issue is settled by Feynman's e parameter, devised to lift infinitesimally the 
degeneracies of the action. It is easy to see that it gives a weak variation of the action along gauge orbits by breaking 
gauge invariance. As long as our truncation of the evolution equation, the functional C[A\ is explicitly gauge invariant 
e plays the role of an infinitesimally external magnetic field in the Ising model, ie helps the breakdown of the gauge 
symmetry only if it is really broken in the true vacuum. 



IV. WHAT HAS BEEN ACHIEVED 



The traditional implementation of the RG idea has proven to be essential in Statistical and High Energy Physics, 
starting with the understanding of critical phenomena [102], finite size scaling [103], e-expansion [104], dynamical 
processes with different time scales [105], continuing with partial resummation of the perturbation expansion [4], 
parametrizing the scale dependence at high energies [106] and ending with the construction of effective theories in 
Particle Physics [107-109]. This is an inexcusable short and incomplete list, its role is to demonstrate variability 
and the importance of the method only. Our main concern here was the functional form of the RG method and 
its generalizations. From this point of view one may distinguish conceptual and more technical achievements as the 
power of the functional formalism is more exploited. 

We should consider the RG method as 'meta-theory', or in more practical terms as a language, as Feynman graphs 
are used in particle physics. But this language is not bound by small parameters and can provide us a general, 
non-perturbative approximation method beyond the semiclassical expansion and numerical simulations. 

The path integration was first viewed as a powerful book-keeping device for perturbation expansion and the truly 
non-perturbative application came later, after having gained some experience with the formalism. This is similar 
to the development of the functional RG method where the RG equation refers to the generator function of the 
effective coupling constants and provides us with a simple procedure to keep track of the Feynman graphs. The 
steps in functional calculus are more cumbersome because they deal with generator function(al)s, with infinitely 
many coupling constants. The ultimate goal is to go beyond this level and to use this formalism in a genuinely 
non-perturbative manner. 

The semiclassical expansion is our first step in this direction. It will be important to check how the RG continues 
to be applicable where the saddle point expansion ceases to be reliable. 

Another field the RG scheme might be compared with is lattice regulated field theory. Both are general purpose 
tools to deal with non-perturbative systems. The bottle-neck of the numerical simulations on the lattice is the need 
to send the U.V. and the IR cut-offs sufficiently far from each others and the restriction to Euclidean space-time. The 
RG strategy is set up in continuous, Minkowski space time and there is no particular problem with keeping the U.V. 
and IR cut-offs far from each other. But the drawback of the RG strategy is that it is rather lengthy to extend the 
space of (effective) action functionals used in the computation. Since the limitations of the two methods are quite 
different they might be used in a complementary manner. 

The functional formalism is promising because of the possibility of following the mixing of a much larger number 
of operators as in the traditional strategy. This feature gives the hope of extending the applicability of the method 
from a single scaling regime to the whole range of scales covered by the theory. Such an extension may provide us 
valuable information about the competition of interactions in realistic theories. 
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